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We numerically simulate ensembles of SU(2) Yang-Mills dimeron solutions with a statistical weight 
determined by the classical action and perform a comprehensive analysis of their properties as a 
function of the bare coupling. In particular, we examine the extent to which these ensembles and 
their classical gauge interactions capture topological and confinement properties of the Yang-Mills 
vacuum. This also allows us to put the classic picture of meron-induced quark confinement, with 
the confinement-deconfinement transition triggered by dimeron dissociation, to stringent tests. In 
the first part of our analysis we study spacial, topological-charge and color correlations at the level 

■ of both the dimerous and their meron constituents. At small to moderate couplings, the dependence 
of the interactions between the dimerous on their relative color orientations is found to generate 

. a strong attraction (repulsion) between nearest neighbors of opposite (equal) topological charge. 

■ Hence the emerging short- to mid-range order in the gauge-field configurations screens topological 
_ I charges. With increasing coupling this order weakens rapidly, however, in part because the dimerous 

, gradually dissociate into their less localized meron constituents. Monitoring confinement properties 

by evaluating Wilson-loop expectation values, we find the growing disorder due to the long-range tails 
of these progressively liberated merons to generate a finite and (with the coupling) increasing string 

■ tension. The short-distance behavior of the static quark-antiquark potential, on the other hand, 

■ is dominated by small, "instanton-like" dimerous. String tension, action density and topological 
susceptibility of the dimeron ensembles in the physical coupling region turn out to be of the order 
of standard values. Hence the above results demonstrate without reliance on weak-coupling or low- 
density approximations that the dissociating dimeron component in the Yang-Mills vacuum can 
indeed produce a meron-populated confining phase. The density of coexisting, hardly dissociated 
and thus instanton-like dimerous seems to remain large enough, on the other hand, to reproduce 

O . much of the additional phenomenology successfully accounted for by non-confining instanton vacuum 

■ models. Hence dimeron ensembles should provide an efficient basis for a rather complete description 
of the Yang-Mills vacuum. 
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I. INTRODUCTION 



The notorious complexity of the Yang-Mills vacuum manifests itself in emergent infrared phenomena which, despite 
vigorous theoretical efforts over the last decades, remain only fragmentarily understood. Quark confinement [1,] was 
early recognized as a paradigm for these infrared complexities, both because of its unprecedented nature and of 
■ its almost universal impact on the hadronic world. The long struggle to explain the nonperturbative confinement 
^ [ mechanism has led to the development and often ongoing refinement of a wide variety of theoretical approaches [H , 
ranging from various vacuum models to ab-initio lattice simulations Q. The latter, in particular, have provided strong 
^ _ and unbiased evidence for the main signatures of quark confinement, i.e. the linear growth (and saturation) of the 
k> . static quark-antiquark potential at large distances and the associated color flux-tube formation [1| (and breaking [5] ) . 
; I ' Beyond establishing confinement and its characteristics, intense efforts were devoted to pinning down the underlying 
, dynamical mechanism and to identify the infrared degrees of freedom best suited to describe it in humanly fathomable 
terms. Vacuum model analyses have played a prominent and often pioneering role in this endeavor. The majority 
of them is based on ensembles of specific "constituent" fields, i.e. at least partly localized gauge fields (with their 
prospective monopole loop and center vortex content) which are supposed to populate and disorder the Yang-Mills 
vacuum. Potential candidates for these building blocks can be classified according to the increasing number of 
spacetime dimensions in which they are localized, and thus equivalently to their increasing efficiency in disordering 
the vacuum. The minimally, i.e. in just two dimensions localized candidates are center vorticeslg, MK An intermediate 
position occupy both (gauge-projected) Abelian '§', '9'] and non-Abelian monopoles or dyons [ifl, lll| (possibly in the 
guise of BPS constituents of KvBLL calorons ^3|], see below) which are localized in the three spacial dimensions. 

Among the maximally, i.e. in all four (Euclidean) spacetime dimensions localized candidates for the buil ding 
blocks of confining field confi gura tions, finally, are regular-gauge instantons calorons with nontrivial [ij, [Tsl - 
[l7i] holonomy and merons [l^. Tisl. [l9| . All of these "pseudoparticles" carry a nontrivial topological charge, mediate 
vacuum tunneling events and solve the classical Yang-Mills equation (ssj . The latter property, in particular, has raised 
hopes that at least one of these solutions may provide a basis for semiclassical and therefore potentially analytical 
treatments of confinement. (This proved indeed possible in lower-dimensional models where instantons or vortices were 
shown to generate confinement at weak coupling [2l|,[22|.) In classically scale-invariant, four-dimensional Yang-Mills 
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theory, however, where one is inevitably faced with large couplings at large distances, the justification for semiclassical 
approximations depends crucially on the characteristic distance scales involved. A successful semiclassical description 
of confinement, in particular, would require all physically relevant confinement effects (as, for example, the quark 
binding inside light hadrons) to take place over distances which remain small enough to avoid the strong-coupling 
regime. 

In order to overcome these potentially stymying limitations, one may resort to simulating fully interacting pseu- 
doparticle ensembles numerically. The first and still most extensive such calculations were based on rather dilute 
superpositions of instantons and anti-instantons in singular gauge [2^ . Such superpositions maintain part of the 
instantons' semiclassical nature, are known to generate a variety of important physical effects including spontaneous 
chiral symmetry breaking, and describe a large amount of successful vacuum and hadron phenomenology [23]. How- 
ever, they fail to generate a confining, i.e. linearly growing potential between static quarks |24l426| . By now, numerical 
simulations were performed for superpositions of all the other above-mentioned pseudoparticle candidates as well. A 
central result was that each of them proved capable of generating a finite string tension. 

The discovery of new finite-temperature instantons [l3l | has given renewed impetus to semiclassically inspired 
confinement models. These generalized calorons contain Nc self-dual monopole ("dyon") constituents without classical 
interactions and carry a non-trivial holonomy as an explicit degree of freedom. The ability of the calorons to separate 
into Nc quasi-independent constituents (without cost in action) in the confinement phase has led to simulations 
[l7| where the separation appears as an additional (internal) degree of freedom compared to the "old-fashioned" 
instanton simulations. Encouraging results have been obtained. Later, an analytic solution for the caloron-dyon gas 
was proposed at the price of a non-positive weight of an overwhelming multitude of configurations with respect 
to the moduli metric [27| . An analytic as well as numerical treatment of the purely Abelian non-interacting dyon gas 
was shown to provide a confining force (28| . 

In addition, equilibrated lattice gauge-field configurations were successfully searched for the presence of pseudopar- 
ticles (with the notable exception of merons) in the Yang-Mills vacuum. Such searches are typically performed by 
smoothing techniques designed to filter out infrared gluon fields and in particular solutions of the classical Yang-Mills 
equation. For examples in the most extensively studied case of instantons and their size distribution see Refs. |29j , 
and for calorons with nontrivial holonomy Refs. jlGj . 

In the present paper, building upon the above developments, we are going to introduce and study a new type of 
pseudoparticle ensemble. The latter is based on superpositions of dimerous and anti-dimerons (in singular gauge), i.e. 
on those four-dimensionally localized solutions of the classical SU(2) Yang-Mills equation which contain two meron 
centers [l8[. The relatively large dimeron solution family includes as one extreme members which are contracted to 
(singular-gauge) instantons and as another members which are dissociated into two far separated merons. Moreover, 
it contains continuous sets of dimeron fields which interpolate between these extremes. Hence our configuration space 
includes many of the previously studied singular-gauge instanton and single-meron configurations and is ideally suited 
to study transitions between instanton and meron ensembles. This is a crucial advantage because such a transition 
process was conjectured to take place between the instanton, dimeron and meron components of the Yang-Mills 
vacuum as a function of the bare gauge coupling. Indeed, this transition is a central ingredient of the meron-induced 
confinement scenario which was put forward on the basis of qualitative arguments by Callan, Dashen and Gross 
(CDG) [2 4| . | 3 0l| (and further explored in several promising directions, including e.g. a two-phase picture of hadron 
structure [3l| which recieved some support from lattice simulations [32|). (For more recent work on merons see Refs. 
[33| . An alternative meron confinement mechanism was suggested in Refs. [3^.) 

The CDG confinement scenario was inspired by the already mentioned weak-coupling confinement mechanism in 
2-1-1 dimensional Yang-Mills- Higgs theory [2l|,[22. In contrast to this classic paradigm, however, even the ingenious use 
of semiclassical arguments turned out to be of limited use in 3-1-1 dimensional Yang-Mills theory. In fact, the resulting 
estimates remained qualitative at best |66| . Further progress was achieved only recently, when a numerical study of 
superpositions of merons and antimerons (l4l. [l9| provided the first convincing support for a central assumption of the 
CDG scenario, namely that the meron component of the Yang-Mills vacuum may indeed confine, i.e. generate a static 
quark potential which rises linearly at large interquark separations. The suggested origin of the meron population 
in the vacuum from dimeron dissociation remained obscure, on the other hand, since the field content was from the 
outset restricted to a purely meronic phase. 

To search for the development of such a meron-dominated phase and to analyze the underlying mechanism in 
dimeron-antidimeron ensembles evolving under the Yang-Mills dynamics will therefore be important points on our 
agenda. For this purpose, and in contrast to previous simulations of pseudoparticle ensembles, we will monitor 
pertinent vacuum properties in a wide range of gauge coupling values. In particular, we will look for signatures of the 
envisioned dimeron dissociation process. Similar to the gradual break-up of molecules into their atomic constituents 
with increasing temperature, this dissociation is suggested to be driven by a competition between the decreasing 
attraction among the dimeron's regularized meron centers and the entropy j67| gain from their increasingly distant 
positions. Guided by the behavior of an ideal dimeron gas, one may indeed roughly estimate the inter-meron attraction 
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to decrease relative to the entropy with growing couphng. In view of the strong interactions anticipated especiahy 
between the meron centers of different dimerons in significantly dissociated ensembles, however, such estimates are 
unreliable. Moreover, in our dynamical context the coupling-dependent competition between energy and entropy is 
particularly subtle since both attraction and entropy are expected to depend only logarithmically on the intermeron 
separation [68i] . This further impedes even qualitative analytical estimates of dimeron ensemble properties. 

Our fully interacting ensembles, on the other hand, are well suited to tackle and clarify these issues. Since they 
take all classical gauge interactions - including those which may be strong, long-range or many-body - between 
the dimerons into account, we do not have to rely on the semiclassical or any other weak-coupling or low-density 
approximation. Nevertheless, it will be of heuristic value and of help for future, possibly more analytic treatments 
to understand which aspects of the dimeron-ensemble dynamics may be described semiclassically. In fact, one of our 
reasons for transforming the individual dimerous into a singular gauge was to provide additional insight into this issue 
by more strongly localizing them. This allows the dimerons and their mcron centers to retain their identities and their 
classical shapes to a larger extent, even in the presence of quantum fluctuations, and thus facilitates (semi-)classical 
behavior. At several points during the course of our investigation we will therefore check for indications of such 
behavior. 

Besides studying the dimeron dissociation process and its dynamics in detail, we will also survey other structural and 
topological properties of the ensemble configurations. In particular, we will examine distributions and correlations of 
the topological charge carriers and search for ordering tendencies with respect to both dimeron and meron centers. (A 
side benefit for investigating topology distributions is that dimerons, in contrast to e.g. center vortices and monopoles 
[ssj . carry the topological charge of the Yang-Mills gauge group directly.) This will result in a more detailed 
understanding of the restructuring processes which accompany the transition to the strong-coupling regime. We will 
pay particular attention to the behavior of the topological susceptibility and of confinement properties, which we 
monitor by evaluating Wilson-loop expectation values and the associated string tension, as a function of the bare 
coupling. The quantitative picture emerging from these investigations will reveal, in particular, at which stage of 
dissociation (as quantified e.g. by the average separation of the dimerons' meron partners) dimeron ensembles can 
best describe the confining phase of the Yang-Mills vacuum. 

In the next section we start by acquainting the reader with pertinent properties of the dimeron solutions, transform 
them into a singular gauge and discuss the regularization of their singularities. We then set up our ensemble field 
configurations and write down their partition function. Section Hill presents an outline of our simulation strategy and, 
in particular, of the measures taken to control systematic finite-size and discretization errors. The next and central 
Sec. IIVI contains our results and their discussion. We start with an extensive statistical analysis of the spacetime, 
topological-charge and color structure of the ensemble configurations as well as of the crucial dimeron dissociation 
process. These investigations provide a broad spectrum of insights into the behavior of the dimeron ensemble and 
its phase structure as a function of the gauge coupling. They furthermore relate this behavior to the changes in 
the average properties of the individual dimerons. We then evaluate the topological susceptibility and the static 
quark potential (based on the calculation of Wilson-loop expectation values) and analyze the results in relation to 
the coupling-dependent restructuring of the ensemble configurations. We further discuss scale-setting issues, evaluate 
dimensionless observable ratios, rewrite our results in physical units and confront them with those of other approaches. 
Section [Vj finally, contains a summary of our main findings and presents our conclusions. 

II. SINGULAR-GAUGE DIMERON CONFIGURATIONS AND THEIR DYNAMICS 

In the following subsections we construct the (continuum) dimeron field-configuration space on which our study 
will be based. Along the way, we will motivate our expectation that such dimeron configurations approximate crucial 
features of the gauge-field population in the Yang-Mills vacuum, and in particular that they may provide a decent 
description of the confinement-deconfinement transition with decreasing gauge coupling. We start by deriving the 
building blocks of our fields, i.e. the individual dimeron solutions of the Yang-Mills equation in singular gauge, and 
discuss their collective coordinates. We then superpose these solutions to define our gauge-field configuration space 
and write down the partition function which specifies its dynamics. Finally, we discuss the (both gauge-dependent 
and -independent) singularities of our configurations and introduce suitable regularization procedures to prepare for 
their numerical treatment in Sec. IIIIl 

A. Yang-Mills dimerons in singular gauge 

Dimerons (or meron pairs) are those classical solutions of the Euclidean Yang-Mills equation which contain two 
meron centers where their only singularities are located [18^. In regular gauge, with the two merons located symmet- 
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rically at the distances ±a from the center at xq, the SU(2) dimeron solution family reads [1| 



4°''-) ix;{xo,a,u})^ 



{x — xo + af' (x — xo — of' 



(1) 



where cr^^ := ij^^^cFa/'^ with the Pauli matrices (Ja and the 't Hooft symbols t]^^^. (The latter are defined as 
= £api/, — ^Va4fi = ^af^ for /i, = 1, 2, 3 and 7]^44 = [l^.) The unitary matrices u are global SU(2) color 
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rotations. The dimerous ([T]) carry the same topological charge as instantons, i.e. Q — 1. The anti-dimeron A 
with opposite topological charge Q = — 1 is obtained from Eq. ([!]) by replacing r]^^^ — ^ f}^^^ := 
contrast to instantons, however, dimerous are not selfdual. 

The solution class ((Ij depends on eleven real and continuous "collective coordinates" or "moduli" whose values 
uniquely and completely specify each member. Eight of these parameters, xq,^ and a^, determine the spacetime 
position and orientation of the dimeron while the remaining three determine the SU(2) group elements u [70|. The 
collective coordinates correspond to those eleven independent combinations of the classical and continuous Yang-Mills 
symmetries, i.e. of Euclidean spacetime translations and rotations, conformal transformations and global SU(2) color 
rotations, which transform a representative dimeron into a gauge-inequivalent solution j40| . 

It is instructive to consider several limits of the solution family ([IJ . For a — )■ the meron centers coalesce and the 
dimeron turns into a pointlike, regular-gauge instanton. (A finite-size instanton is reached after suitably regularizing 
the meron singularities, cf. Sec. Ill CI ) For \a\ — >■ oo, on the other hand, the dimeron breaks up into two merons. 
Although these merons maintain the fixed relative color orientation of Eq. ([1}, changing it by hand will yield 
increasingly action-degenerate two-meron solutions when the inter-meron separation \a\ becomes large. This is because 
the color-dependent attraction which locks the meron constituents of Eq. ([I} into their rigid color orientation decreases 
with increasing \a\. Hence both instantons and merons are contained in the solution class ([1]) as limiting cases. The 
fields (HI therefore provide on-shell interpolations between instantons and isolated meron pairs, i.e. continuous paths 
in Yang-Mills solution space which connect these particular gauge fields. In our context these paths are of particular 
interest since they provide preferred doorways along which the dimerous may dissociate into merons. As already 
mentioned, such dissociation processes are conjectured to drive the deconfinement-confinement transition. 

The solutions (H]) have the same \x\ — oo behavior as an instanton in regular gauge [H, HJIi^, i.e. they contain a 
long-distance tail 
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~ 1/ I a; I which implies an exceptionally weak localization. In multi-dimeron configurations these tails generate strong 
overlap interactions between the individual dimerous. Building dimeron ensembles by superposing solutions of the 
type ([1]) would thus lead to very strongly correlated systems, in some respects similar to the regular-gauge instanton 
ensembles studied in Ref. [l^l. In the present work we will take a different route, however. In contrast to merons, 
dimerous - like instantons - may be transformed to singular gauges in which they are more strongly localized, with 
their long-range tails decaying as A superposition of singular-gauge dimerons thus improves the vacuum 

description at short distances, compared to superpositions of either regular-gauge instantons or single merons. As 
long as the characteristic dimeron size is small compared to the average inter-dimeron separation, furthermore, it will 
provide much better approximations to classical Yang- Mills solutions. This is essentially because at moderate pseu- 
doparticle densities the contribution of the overlap regions to the action is much smaller than in regular-gauge dimeron 
superpositions. Hence singular-gauge multi-dimeron configurations provide a privileged testing ground especially for 
the semiclassical features of the CDG confinement mechanism [23. [30l[3l|. 

In addition, dimeron configurations in singular gauge may also describe less classical or even fully quantum- 
mechanical aspects of confinement. This holds, in particular, for the vacuum disordering mechanism envisioned 
by CDG. The latter relies on the long-distance tails of the individual merons which come into play when the separa- 
tion \a\ between the meron centers of the dimerons becomes large. As discussed above, the partner merons are then 
practically independent and regain their long-range tail 1/ |a:| (at least in isolation). The ensuing long-distance 
color correlations among these merons were found in Refs. |14l . [l9j to be far from semi-classical and to generate 
confinement. Of course, a complete dimeron dissociation into isolated and thus far separated merons is impossible in 
the finite spacetime volumes in which numerical simulations are feasible. (The same holds for the solution in which 
one of the merons delocalizes on the three-sphere at spacetime infinity ^ll,[23|-) However, such isolated merons would 
anyhow be unphysical (since they carry infinite action) and incompatible with a finite dimeron density. Instead, one 
expects that beyond some typical separation the two meron centers of a dimeron experience stronger interactions with 
their individual field environment than with their partner. The rigid link between their color orientations can then 
be broken, i.e. the partner merons can become effectively independent of each other. 
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In order to construct the multi-dimeron configurations motivated above, we first transform tfie dimeron solution 
family ^ into a specific singular gauge, 

4°'^) = f/t (x) A^^'^^Uix) + *C/t (x) d^U (x) , (3) 

where U is the large, i.e. topologically active gauge group element 

U{x)=tx^al+^ (4) 

(tJ^'^-' = (a, —i) , x^ = Xf^/ y/x^x^) which has a singularity at the origin of the coordinate system. The result can be 
cast into the form 

4°'^-' = /i (x, a) x^a^^ + /a (x, a) a^X^,^ {x) (5) 
{^^lv '■— Vanu^o./'^)- Choosing for simplicity xq — Q and u — 1, the two scalar functions /i_2 become 

/i {x, a) := ^ - - — - - — 5— 2 , (6) 

x^ [x + a) {x - a) 

h{x,a) ,^2 ~ . ^ .2 - (7) 

(a; + a) [x — a) 

The antisymmetric and anti-selfdual field X^^ contains important parts of the SU(2) and spacetime tensor structures 
(and their mutual couplings). Explicitly, 

Xf_,^ (x) = Va^ii^Xa (x) with Xa {x) = Q ~ ^r^r^ CTa + XaXsCTs - X4XsSascCrc (8) 

(latin indices are spacial). After regularization of the short-distance singularity (cf. Sec. Ill Cp the a limit turns 
the dimeron (O into a singular-gauge instanton with /i (a;, a = 0) = 2p^/ [x^ (a;'^ -I- p^)] and /2 (a;, a = 0) = 0. In 
the following, we will often refer with the term "dimeron" to both dimerous and anti-dimerons in singular gauge. 
Occasionally, we will use the term "pseudoparticle" for the same purpose. 
The leading asymptotic behavior of the singular-gauge dimerous ([5|) is 



i.e. the long-range tail © has as intended disappeared and the overlap between neighboring dimerons is strongly 
reduced. (This is in contrast to the meron ensembles of Refs. [13, [13 whose constituents exist only in regular gauge.) 
In addition to the meron-center singularities at x = ±a, the solution ([SJ inherits another singularity at a; = from 
the gauge transformation Q . Hence the impact of the latter will disappear when forming (topologically insensitive) 
gauge- invariant quantities from Eq. ([S]). The regularization of these singularities will be discussed in Sec. Ill CI 

By construction, the singular-gauge dimerons ([5]) provide only a subset of the full solution family. The complete 
eleven-parameter family can be recovered by translating the solutions ([SJ to xq and by gauge-rotating them with a 
constant matrix u G SU(2), i.e. 

{x; {xo, a, u}) = ^t^P.s) _^^)u = R^b {u) Af^'^^" (x - xq) y (10) 

where the orthogonal matrices Rab {u) = tr {craitVbw} /2 are obtained from the condition UbU = Rab {u) da- The 
group elements u depend on three real and continuous parameters which form (local) coordinates on the SU(2) group 
manifold . In the following it will be convenient to use the quaternion representation 

U (c) = Co + iCaCTa (11) 

which embeds into with Euclidean coordinates co, C3 by imposing the unit three-sphere constraint X)a=o '^a ~ 
1. In these coordinates R takes the form 

Rab {u (c)) = (co - CcCc) Sab + "^.CaCb - 2eabcCoCc. (12) 

For our discussion below it will be useful to keep in mind that dimerons have three more continuous and noncompact 
collective coordinates than instantons (but three less than a pair of independent merons, due to the locking of the 
inter-meron color orientation). The three additional coordinates arise from the more complex structure of extended 
(i.e. \a\ ^ 0) dimerons which requires more degrees of freedom to locate and orient them in spacetime. This results in 
a substantially larger "position entropy" [2^ |43| which is instrumental in counterbalancing the (regularized) dimeron 
action which grows logarithmically with \a\. Indeed, the larger entropy is a necessary (but not sufficient) requirement 
for dimerons to dissociate with increasing coupling and to finally split up into their meron partners. 
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B. Field configurations and partition function of the dimeron ensembles 



As motivated above, we intend to study a model which drastically reduces the field content of Yang-Mills theory 
(as integrated over in amplitudes and the partition function) to superpositions of Nd dimerous and Nq ~ Njj 
antidimerons in singular gauge [71|, i.e. 

No No 

A, (x, {C,,„C,- J) - ^4°'^) (x; {Cj4) + 5^4°'^) (x; {C%}) (13) 
/ I 

with TV = Njj + Nf). Each term in this sum is uniquely characterized by the set {Cj} = {xQj,aj,cj} of collective 
coordinates of the corresponding (anti-) dimeron. We recall that the configurations ([T3l) differ distinctly from those 
obtained by transforming a regular-gauge dimeron superposition into (any) singular gauge. This is because in Eq. (|13p 
the gauge of each dimeron is chosen relative to its individual position. One may wonder, incidentally, whether n-meron 
solutions with Q > 1 and their anti-solutions should be added to the superposition ansatz (jlSp . As in the case of 
instantons, however, the Bogomoln'yi-type bound S > [Sn"^ /g^) \Q\ on the action [4l| implies that such multi-(anti)- 
meron contributions to the partition function are exponentially suppressed relative to the dimeron contributions 
(t^ . Since the entropy increases only logarithmically with n and thus cannot compensate this suppression, such 
multi-meron contributions may be safely neglected. 

Nevertheless, the ansatz (jl3p should be regarded as a rather minimal choice. It is mainly geared towards a trans- 
parent study of the proposed (instanton and) dimeron dissociation mechanism and its role in the deconfinement- 
confinement transition [23 . [sol [sij . Hence there are several natural directions in which Eq. (jlSp may be extended in 
future studies to provide a more complete description of the Yang-Mills vacuum physics. An example would be to 
add the meron-antimeron pair solutions [l^ of the Yang-Mills equations, again individually transformed into singular 
gauge. This would maintain the approximately semiclassical nature of the configurations at small \a\ and yield a 
richer dynamics. However, it would probably also lead to a less transparent interpretation of the results, and we do 
not expect the topologically trivial meron-antimeron pair configurations to provide qualitatively new insights into the 
transition behavior. Indeed, their limits for \a\ — >■ 0, a pure-gauge field of zero action, and for \a\ — )■ cxd, a meron 
and an independent antimeron, indicate that meron-antimeron pairs do not generate new pathways for the transition. 
More promising improvement options would include generalizations of the dimeron superposition ansatz (1131) which 
allow for a complete transition into a meron ensemble (e.g. by releasing the rigid color locking between the dimerous' 
meron partners beyond a suitable intermeron separation lajj^^jj^) or the admixture of a pure instanton component 
[13, [llj with a realistic size distribution [iS] . 

As discussed above, we view the dimeron configurations ([T3l) as a pertinent subset of the SU(2) gauge fields governed 
by the Yang-Mills dynamics. Hence we define the partition function of our dimeron model as 

.NoMo 

Z= n dCidCre^p{-S [A^{{Ci^,,CrJ)]} (14) 

where S is the Euclidean Yang-Mills action 

5 [^] = ^ y" d^xtr {F^.F^,} J d^xs [x) (15) 

and Ff^u is the gauge field strength 

{x) - d^,A, {x) ~ d,A^ (x) - i [A^ (x) , A, (x)] . (16) 

As already indicated, we expect the gauge interactions among the dimerons in Eq. (|14p to play an important role in 
generating confining long-range correlations for large \a\ and, in particular, to provide a disordering mechanism for 
the vacuum. This is in contrast to the situation in non-confining singular-gauge instanton ensembles where for most 
amplitudes a random orientation of the instantons (i.e. the neglect of the interaction term exp (—5) in Eq. (|14p ) 
yields a reasonable approximation to the vacuum physics [23j . 

The integration over the collective coordinates will be performed with the measure 

dC = d^xod'^adfi {c) (17) 

{dfi (c) is the SU(2) Haar measure) for the dimerons, and analogously for the anti-dimerons. This type of measure 
if familiar from instanton vacuum models (23| . An improved alternative may include Jacobians which arise from the 
transformation of the linear gauge-field measure dA into the collective-coordinate basis and thereby implement the 
full moduli-space metric. Finally, an additional a dependence of the measure can emerge from the trace anomaly 

MM- 
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C. Regularization of the gauge- field singularities 

The singular-gauge dimeron solutions ([5]) contain two types of singularities. The first are those located at the 
constituent meron centers x = Xq zL a which persist in any gauge. In addition, there is a gauge singularity inherited 
from the topologically large and thus necessarily singular gauge transformation 1^, chosen to sit at the origin. 
In isolated dimerous, such gauge singularities could be gauged away and would therefore not affect (topologically 
insensitive) observables. This ceases to be the case for the superpositions p^ . however, where the location of the 
gauge singularities varies with the positions of the pseudoparticles. 

Since both types of singularities would impede our numerical simulations, they have to be regularized in a physically 
acceptable manner. (Recall that the individual dimeron fields Q and ([5]) do not solve the Yang-Mills equation at 
their singularities, so that a suitably localized regulator will only minimally affect the semiclassical properties of the 
configurations p3ll .) A natural way to regularize the singularities at the meron centers is to add the square of a 
"size" parameter p to the denominators {x — xq i: of' in Eqs. ([6]) and ([7]). Such a regulator imitates the way in 
which the size parameter enters the instanton solutions. Hence it may arise from scale-symmetry breaking quantum 
fluctuations which are expected to smear the classical singularities. A finite p furthermore acts as a UV cutoff since 
it limits the gradients of the configurations (IT3t . Of course, the regularized field configurations cease to be exact 
solutions of the classical Yang-Mills equation. Hence the value of p should be chosen large enough to avoid sizeable 
discretization errors (cf. Sec. IIII Pp but also small enough to avoid unnessary deformations of the dimeron solutions 
from the semiclassical saddle points. We have found p^ = 0.3 to be a reasonable compromise between both criteria 
and will use this value in our numerical simulations, if not stated otherwise. 

The remaining singularities are, at least for isolated dimerous, gauge artefacts. They remain unlikely to have a phys- 
ical impact in our dimeron superpositions, too, but they may nevertheless cause problems in numerical approximations 
(even for single dimerous), e.g. when subtle cancellations in gauge-invariant quantities are upset by discretization 
errors. A finite-size regulator as above would hardly help to avoid such problems since it would just spread out 
the unphysical action-density peaks at the singularities. Hence we treat these gauge-dependent singularities in a 
pragmatic way, namely by interpolating the gauge field in 4c? balls of radius e around the singularities with the field 
value at a specific point on the surface of the ball. In practice, it turns out that with — 0.3 these highly localized 
singularities are smoothly regularized by taking e as small as £ = 5 x 10""'^'^. (The values of dimensionful quantities, 
such as p and e above, are given in "numerical units" originating from the discretization grid to be introduced in Sec. 
IIII Al As a convenient length unit we have chosen a = 0.1, i.e. a tenth of the distance between nearest grid points. 
(Of course, this a should not be confused with the four-vector a which parametrizes the dimeron solutions.) When 
transforming our results into physical units starting from Sec. IIV C[ we will denote quantities in these numerical units 
by a hat above their symbols, cf. Sec. IIV CI ) In our simulations it typically took several hundred sweeps through 
configurations of = 487 pseudoparticles (cf. Sec. IIIIB|) before such a singularity was first encountered. (Even with 
a dramatically reduced floating-point variable length of 16 bits, incidentally, we have observed no overflows (which 
unregularized gauge singularities would generate) when evaluating the action density of a regularized two-dimeron 
configuration.) 

III. SIMULATION DETAILS 

In the following subsections we summarize how we have generated the dimeron configuration ensembles on which 
all our subsequent calculations will be based. We further discuss predominant sources of systematic errors in these 
ensembles. 



A. Sampling volume and resolution 

Owing to the translational invariance of the Yang- Mills dynamics ([T5|) . ensemble averages over the pseudoparticle 
configurations (jl3p in infinite spacetime will likewise be translationally invariant. In our simulations we are restricted 
to a bounded volume of numerically manageable size, on the other hand, and thus have to keep boundary artefacts 
in calculated amplitudes under control, ideally within the size of the statistical uncertainties. Our initial step in this 
direction was to adopt singular-gauge dimerous with their more rapidly decaying long-range tails as the constituents 
of our gauge-field configurations (IT^ . A sufficient suppression of boundary artefacts turned out to require additional 
measures, however, since the overlap among the still rather moderately localized dimerous generate important (and 
typically repulsive) interactions over ranges beyond the average nearest-neighbor distance. In fact, these interactions 
are expected to play an important role in the confinement mechanism and should thus be distorted as little as possible 
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by the boundary. Simulation costs, on the other hand, with their by far largest part due to integrating the action 
density on a spacetime grid, should be kept minimal. 

In order to approximately meet these conflicting goals, we will adopt a multi-layered description of the boundary 
with decreasing grid-point density towards the outer layers, as sketched in Fig. [TJ The innermost "core volume" is a 
rectangular spacetime box in which we intend to evaluate amplitudes and thus require the physics to best approximate 
the infinite-volume limit. In order to allow for an efficient evaluation of Wilson loops with elongated time directions 
in Sec. IIVFI this volume and its grid are chosen to be asymmetric. One dimension, singled out as the Euclidean time 
direction, contains 40 grid points while each spacial dimension contains 25. This turns out to provide a sufficiently 
high resolution for the amplitudes to be calculated below (cf. Sec. IIIID|) . 

The core volume is encompassed by an "ensemble volume" which contains the (anti-) meron centers of all (anti-) 
dimerous. Its surface is implemented numerically by rejecting Metropolis updates (cf. Sec. IIII Bl) during which a 
meron center would leave this volume. The extent of the ensemble volume in a given direction will be taken 15% 
larger than that of the core in the same direction. In the part of the ensemble box which surrounds the core we reduce 
the grid-point density to (2/3) ~ 20% of that in the core volume. This turns out to yield a still adequate resolution 
while substantially reducing the computational cost of evaluating the action. 

A satisfactory suppression of field distortions inside the core volume turns out to require an additional precautionary 
measure, however. It consists in correcting for a particularly prominent boundary effect, namely the artificial attraction 
of the pseudoparticles to the surface of the ensemble volume. The latter arises because a substantial part of the action 
density of dimerous near the boundary is located outside the ensemble volume and thus not accounted for while 
the compensating tails of outside dimerous reaching into the ensemble volume are neglected. We approximately 
remove this artefact by surrounding the ensemble box with another, "covering" volume which extends beyond the 
ensemble box by 10% of the core size in each direction. To keep the additional computational costs under control, we 
decrease the grid-point density of the outermost shell (i.e. the part of the covering volume not shared by the ensemble 
volume) inversely with the distance from the ensemble boundary. Since this shell does not contain the rapidly varying 
fields close to the meron centers, calculating the action in the full covering volume indeed largely prevents the fake 
attraction to the ensemble boundary (cf. Sec. IIVB[) . Alternatively, finite-size effects in pseudoparticle ensembles can 
be efficiently corrected by employing the Ewald summation technique, as has been shown very recently in the simpler 
case of Abelian dyon field ensembles psj . 

B. Monte-Carlo updates with dynamical resolution and step-size adaptation 

We evaluate the discretized functional integral over the pseudoparticle fields, i.e. the multi-dimensional integral 
over their collective coordinates in the partition function ([T^ or any other amplitude, stochastically by Monte-Carlo 
importance sampling. Hence we average over dimeron configurations chosen randomly from a Gibbs distribution with 
Boltzmann factor exp (—5*) where S is the Yang-Mills action (fT5|) . More specifically, we use the Metropolis algorithm 
to generate homogeneous Markov sequences of dimeron configurations that visit fields with larger probabilities more 
often. After reaching equilibrium, the probability of finding a configuration in the ensemble of subsequently generated 
fields follows the Gibbs distribution. 

The initial pseudoparticle configurations for these Markov sequences are obtained by choosing their collective 
coordinates from a uniform random distribution. This procedure has equal access to all members of the configuration 
space and practically always results in configurations far from equilibrium, with actions several orders of magnitude 
above their equilibrium values. We have also tested various more ordered initial arrangements of the dimerous and 
convinced ourselves that those lead within errors to equivalent equilibrium ensembles (see Sec. IIII Pp . 

For the Markov update of the n-th configuration An to its sequel An+i according to the Metropolis rules, we first 
generate some candidate configuration A' by randomly choosing an individual pseudoparticle in An and modifying 
alternately (from one candidate to the next) either its position, i.e. — > + da^t and xq.^i — )■ xq.^i -I- dxo,fj,, or its 
color orientation Ci — > -|- dfi (ci). The increments da, dxo and dc are chosen randomly and uniformly subject to the 
constraint that their length remains limited. More specifically, we demand \dxQ\ € [0, Ma^olmaxl' M*^! ^ M'^lmaxl 
well as \dc\ S [0, Irfcl^j^j^^] [tI] and further restrict jdalmax = I'^-^olmax- '^^^ candidate configuration A' is accepted as 
the new configuration An+i with probability min {1, exp (S [A] — S'[A'])}. If a candidate is rejected, the unchanged 
configuration is taken as An+i (and included in measurements like all others). 

The choice of the maximal modification step sizes |(ia|,„„„ and Idcl ,„„ can be used to optimize both the thermal- 
ization rate and the decorrelation among subsequent field configurations of the generated ensembles. This requires 
a compromise between too small values, which impede the progression through configuration space, and too large 
values which more strongly decorrelate subsequent configurations but cause an unefficiently large update rejection 
rate. Examination of the relation between maximal step sizes and typical acceptance probabilities suggests an optimal 
acceptance rate of about 25% for updates in both position and color space. During the initial thermalization process 
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we will therefore increase (decrease) the maximal step sizes Ic^alj^j^^^ and Idcl^j^j^^ after every 10 consecutive candidate 
field configurations by 10% if the average acceptance rate for these configurations falls below (rises above) 25%. 

This dynamical step-size adaptation procedure accelerates the approach to equilibrium and avoids getting trapped 
into approximate would-be equilibria. Initially, i.e. far from equilibrium and near the high-action random configura- 
tions, the fields can relax in larger steps while closer to equilibrium the action fluctuations become smaller and require 
decreasing step sizes to maintain sufficient acceptance rates. From the time when approximate equilibration sets in, 
however, the step size is kept constant to preserve detailed balance. We further note that increasing g increases the 
acceptance probability since the Yang-Mills action ((T5|) scales as g"^. Hence larger g allow for larger step sizes, accel- 
erate thermalization (in less Markov steps) by decreasing the autocorrelation time and result in equilibrium ensembles 
with larger entropy. 

As a case in point, for ~ 1 (and = 0.3) wc find l^alfnax ^ ^-^^ M'^lmax ^ ^-^ equilibrium, where the 
scale of the fluctuations is set by the competition between action and entropy. (For = 0.2, i.e. for more strongly 
localized meron centers, one instead reads off the smaller values Idalj^^g^^ ^ 0.025 and l^clj^^g^^ ^ 0.2 from Fig. [2]) For 
g2 — 10^^ on the other hand, one finds the indeed substantially larger maximal step sizes Malmax ~ [^"^Imax ^-^ 
(again for ~ 0.3). The decreasing step width during a typical thermalization history is plotted in the uppermost 
row of Fig. [2] for g^ = I and g^ — 25. Note that different step sizes in position and color space are generally required 
to obtain action changes of comparable magnitude. 

The efficiency of the ensemble generation process can be further improved by exploiting the during equilibration 
decreasing ruggedness and action of the dimeron configurations in yet another way, namely by starting on coarser 
grids and increasing the resolution only when a higher accuracy of the action evaluation becomes necessary. Indeed, 
as long as the action density is comparable to or larger than its value around the regularized meron centers, a reduced 
resolution is generally sufficient and can save computer time. During thermalization the overall action decreases 
strongly, however, and the grid must be refined to prevent the then more prominent meron centers from artificially 
reducing their action density by "hiding" between grid points. In practice, we implement this refinement procedure 
by starting the simulations with a (1/3)'' times smaller grid-point density. When the values of suitable quantities 
(e.g. the distance between nearest- neighbor meron centers) begin to saturate, this factor is increased to (2/3) . Only 
after again reaching approximate saturation the full grid-point density of Sec. IIII Al is activated for the final approach 
to equilibrium and the subsequent generation of the thermal ensembles. The approximate saturation plateaus and 
subsequent grid refinements (indicated by dotted vertical lines) are clearly visible in the thermalization histories of 
Fig.H 

In the following we will denote a set of 2A^ consecutive Markov steps as a "sweep" (where N = Njj + Njy, cf. 
Sec. IIIBI) . During such sweeps each pseudoparticle of a configuration is on average once considered for a full update. 
(The factor two arises since the individual Markov steps attempt to modify either a dimeron's position or its color 
orientation, i.e. only one of the two subsets of its degrees of freedom.) After four consecutive sweeps with the 
acceptance rate kept fixed at 1/4, all pseudoparticles in a configuration are therefore on average updated once. In 
order to reduce autocorrelations among ensemble configurations and thus to increase the statistical independence of 
successive measurements, we will only select the configurations generated by every fifth sweep (after approximate 
equilibration of the Markov chain) as ensemble members and employ a binning procedure to calculate amplitudes and 
their errors (cf. Sec. IIII E[) . 

C. Approach to equilibrium 

As indicated in the previous section, a reliable calculation of vacuum expectation values as ensemble averages 
requires a sufficient thermalization of the Markov sequences which generate the ensemble configurations. Our criteria 
for when their distribution approximates the equilibrium distribution closely enough are based on the measurement 
of several observables and correlations to be described below. 

First insights into the configurations' thermalization properties can be gained by following the Markov evolution 
of the average distance d of a fixed meron center to its nearest neighbor as a function of the number of sweeps. 
Two examples for such histories are depicted in the second row of Fig. [21 In the left panel the coupling is set to 
g = 1, i.e. our smallest value, which results in the slowest equilibration rates we have to deal with in this paper (cf. 
Sec. IIII Bp . (For these measurements we have also selected a smaller meron-size regulator p^ = 0.2 which further 
slows the equilibration process.) A first and general insight to be read off from this figure is how the degree of 
thermalization reached after a given number of sweeps depends on the measured quantity. While the average distance 
to the closest neighboring meron of equal topological charge appears to have equilibrated after about 800 sweeps, the 
analogous distance for oppositely charged merons may not be fully thermalized even after 1200 sweeps. A similarly 
long equilibration history can be observed for the average probability / that the nearest neighbor of a given meron 
center has opposite topological charge, as plotted in the last row of Fig. [51 This behavior seems to indicate that the 
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average distance between two meron centers from the same dimeron equilibrates faster than their average separation 
from the meron centers of neighboring dimerons. For = 25 the typical relaxation times reduce to at least half of 
those at = 1, confirming the expectation that more strongly coupled ensembles thermalize faster. 

For a better understanding of the global thermalization behavior, and especially of the character and uniqueness 
of the reached equilibria, we have also compared ensembles resulting from different initial configurations. More 
specifically, we have prepared a set of initial dimeron arrangements in which the pseudoparticles were regularly 
positioned at maximal average distances. Their topological charges, color orientations and sizes were chosen either 
to approximate different types of action minima or to imitate close-to-equilibrium configurations found in previous 
thermalization runs. 

All ensembles resulting from these different initializations turned out to generate within errors the same averages. 
This provides, for once, evidence for the ergodicity of the underlying Markov process, i.e. for its ability to access 
any possible dimeron configuration in a finite number of steps. More importantly, this result strongly supports 
the conclusion that the thermalization processes indeed have come sufficiently close to equilibrium, i.e. that our 
measurements in Sec. IIVI are performed in almost thermalized ensembles. This conclusion is strengthened by the fact 
that all our runs with non-random initializations were performed &t ~ 1 where equilibration is particularly slow. 

D. Finite-size and discretization errors 

We now turn to the analysis of the systematic errors which arise from finite-resolution and finite-size effects. We 
then describe the steps taken to control them by accordingly refining our calculational strategy. Additional error- 
reduction measures, which apply to the calculation of specific amplitudes only, will be outlined in their corresponding 
sections below. 

Discretization errors arise in our context from the practical necessity to sample the field configurations (1131) on 
spacetime grids of finite resolution. To meet our accuracy requirements, we adapt the grid-point density according to 
the characteristic length and gradient scales in the pseudoparticle configurations under consideration. During thermal- 
ization, where these scales change drastically, we do so dynamically as described in Sec. IIII Bl In the (approximately) 
thermalized ensembles, on the other hand, these scales are essentially fixed by the pseudoparticle density and by the 
size of the regularized meron-center singularities ff^. As described in Sec. Ill CI the regulator is chosen about five 
times larger than the lattice unit a = 0.1, at p ~ 0.55. This proves sufficient to keep the action of an isolated dimeron 
practically independent of its position on the grid. Although configurations with gradients larger than those around 
single meron centers frequently occur among the dimeron superpositions (1131) . their enhanced action renders them 
relatively unimportant when thermalization is achieved. This explains why we did not encounter such configurations 
in our equilibrated ensembles (cf. Sec. IIVBI) . 

The grid resolution's impact on the calculated action values and on the thermalization process can be seen directly 
in the Markov evolution histories of Fig. [2] The two vertical lines indicate the sweep numbers at which the resolution 
of the grid is refined (cf. Sec. IIIIB|) . The full resolution, corresponding to a = 0.1, is reached only after the second 
refinement step. In order to provide a particularly stringent test of discretization errors, we have reduced our regulator 
p for this simulation by a factor of four, to p = 0.14. This size is of the order of the minimal lattice unit and about 
half of the value p ^ 0.3 below which discretization errors become a concern. As a consequence, one observes that the 
evolution of both plotted quantities (i.e. the average meron-center distances (rfnearestMivi) and the probability for the 
nearest neighbor- meron to have opposite Q = ±1/2) on the coarsest grid begins to saturate at a preliminary would-be 
equilibrium when further relaxation is prevented by insufficient action-density resolution. Already after the first grid 
refinement, however, relaxation continues until the deviation from the thermal values reduces to about 15 - 20%. 
The equilibrium values are approached after the second refinement step. For our standard meron size p ~ 0.55 the 
discretization errors will of course be much smaller and should be well under control. This conclusion is confirmed by 
checks on other observables, e.g. when calculating link elements and their concatenations into Wilson loops in Sec. 

EYB 

We now turn to the discussion of finite-size effects. Those are artefacts of the simulation volume's boundary and turn 
out to be more difficult to control than the discretization errors. Our choice of the more strongly localized dimerons 
in singular gauge ([5]) as the constituents of the configurations (jl3p was partly motivated by reducing such boundary 
effects. For the same purpose we designed the multi-layered boundary outlined in Sec. IIII Al Nevertheless, sizeable 
boundary artefacts of different origins, of different g^ dependence and with varying impact on the calculated quantities 
remain to be dealt with. In order to analyze those quantitatively, we start by monitoring violations of translational 
invariance in the ensemble-averaged action density (s (x)) as a function of the distance from the boundary for two 
values of g^. The results are of direct physical and practical importance since (s) is related to the gluon condensate 
and since a reliable evaluation of the action is indispensable for generating trustworthy ensembles. 

The action density s (x) , as defined in Eq. (|15p , can be obtained analytically by inserting the dimeron superposition 
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([T3)) into the field strength ([16]). The ensemble average (s) is then evaluated as outlined in Sec. IIII Bl with ensemble 
configurations taken from every fifth sweep after equilibration. In order to visualize the local homogeneity properties 
of the action density in different sectors of the simulation volume, we plot (s) in Figs. |3] and |4] as a function of the 
distance from the center of the box along two different directions. The latter were chosen to provide insight into the 
anisotropy of the boundary effects, caused in particular by our asymmetric grid with its elongated temporal direction. 
(Instead of calculating the values of s (x) at neighboring, equidistant points throughout these directions, we surround 
them with adjacent hypercubic boxes of side length 0.2 inside which we sample s (x) at 10 randomly chosen points.) 

We start our discussion of the (s) plots with those generated at our smallest coupling value, = 1. In the left 
panel of Fig. [3] we show (s) (black squares, full line) along the time direction, in the right panel along the average 
over all diagonals of the spacial box at the mid point of the time axis. Deviations from a translationally invariant, 
i.e. constant (s) near the center of the box are somewhat smaller along the time direction than along the diagonal. 
This is expected since the former keeps a larger distance from the spacial boundaries. Nonetheless, (s) fluctuates 
considerably even in the temporal direction and even close to the box center. A reasonable fit to a constant plateau 
remains possible inside most error bars up to distances of order 1.0 both along the time and diagonal directions, 
however, as also shown in Fig. [S] Beyond distances of about ten lattice units the missing field contributions from 
outside of the box volume (cf. Sec. IIII A[) begin to reduce the action density substantially. 

As a consequence of these results, we will restrict the volume in which we evaluate amplitudes to subvolumes of the 
core box in which the maximal distances from the center are of order one (if not noted otherwise). To demonstrate 
the dependence of the boundary effects, we further show the two analogous (s) profiles for g^ — 25 in Fig. U) 
Clearly, the boundary effects are strongly reduced and the regions close to the center now show manifest plateaus 
which extend up to distances ~ 1.5 from the center. With further increasing g^ the boundary artefacts become even 
more restricted to the surface region until the growing dimeron dissociation (i.e. the growing intermeron distance 
2 \a\) creates a different type of sensitivity to the boundary, as discussed in Sees. IIII Al and IIVBI 

Since in general both character and strength of the boundary artefacts are amplitude dependent, we have also 
monitored the vacuum expectation value of quadratic "probe" Wilson loops W as a function of their distance from 
the box center. The results will guide us in reliably calculating the expectation values of rectangular Wilson loops 
in Sec. IIVFI Since the largest among these loops play a crucial role in understanding the confinement properties of 
dimeron ensembles, it becomes especially important to understand the boundary's impact on them. We survey it by 
evaluating the expectation values {W) of quadratic Wilson loops (for details see Sec. IIVFI) of side length 0.5 centered 
along the same two lines as the action density (s) above. The loop orientation turns out to have no significant impact 
on the value of (T^) and will thus be averaged over. In Fig. |3] we have included plots of 450 x (W) along the time 
direction in the left panel and along the average over the spacial diagonals in the right panel for g'^ = 1, and in Fig. 
m we display 20 x (W) for g^ — 25. In all four cases, with growing t,r the distance dependence of (W) becomes 
increasingly proportional to that of the average action density (s) . For quadratic Wilson loops of minimal size this is 
expected since such "plaquettes" form the essential part of a discrete approximation to the Yang-Mills action density. 
In our case, the observed proportionality further indicates that the configurations are sufficiently smooth over our 
rather large test loops, probably because their side length is comparable to the regulated meron size p ~ 0.55. 

To summarize, while discretization errors in our dimeron simulations can be reliably controlled, we have observed 
significant violations of translational invariance in both action density and Wilson loops. Although our efforts to 
reduce boundary artefacts prove effective close to the center of the core box, substantial deviations from homogeneity 
set in at distances of order one, in particular for our smallest coupling value g — 1 where they are most pronounced. 
These boundary effects originate from a combination of the strong intermediate-range interactions between dimerous 
and the collective effects emanating from the boundary (cf. Sec. IIII A[) . The lessons learned from the above analysis 
will guide us in keeping boundary artefacts of our results tolerably small, mainly by relying on specifically reduced 
volumes in which to evaluate amplitudes. 

E. Ensemble statistics 

The effective dimeron field theory introduced in Sec. [Ill has three adjustable parameters: the gauge coupling g^, 
the meron-center size p, and the approximately equal numbers N]j ~ Njy of dimerous and antidimerons. We have 
generated all our ensembles with N = Nu + Njy — 487 pseudoparticles, divided into Nu = 243 ~ N/2 dimerous and 
Njy = 244 ~ N/2 antidimerons, in the multi-layered multi-grid volume specified in Sec. IIII Al The meron centers 
are smeared as described in Sec. Ill CI with a common regulator value = 0.3. We have generated two independent 
ensembles for each of the coupling values g^ € {1, 25, 100, 1000, oo} (with decreasing number of members) and a third 
one for g^ = 1. The total number of sweeps per coupling in equilibrium as well as the maximal stepsizes are collected 
in Tab. U 

In order to estimate autocorrelation effects in the ensembles with g^ — 1, where they should be maximal, we 
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2 

9 


# of sweeps 






1 


1843 


0.04 


0.11 


25 


1550 


0.20 


0.40 


100 


1470 


0.36 


0.44 


1000 


554 


0.47 


0.49 


oo 


574 


0.10 


1.00 



TABLE I: Dimeron ensemble characteristics. (For definitions see Sec. IIIIBI ) 

have calculated autocorrelation functions for several quantities of interest, including the average distance between 
nearest-neighbor meron centers of equal and opposite topological charges. For those of opposite topological charge 
the autocorrelations decay fastest, after about 300 - 400 equilibrium sweeps, while this takes up to a few hundered 
sweps longer for equally charged neighbors. The autocorrelation functions for the average distance between nearest- 
neighbor dimerons, on the other hand, show no obvious preference for either equally or oppositely charged neighbors. 
Both distances approximately decorrelate after 500 sweeps. The same holds for the average distance between the 
dimerous' meron partners. 

Based on the above observations, we have adopted the following strategy to reduce autocorrelations. The ensembles 
contain the configurations generated by every fifth Metropolis sweep. The total number of equilibrium sweeps (per 
coupling) is divided into ten bins. (For = 1 a, bin thus contains about 180 sweeps and about 36 ensemble 
configurations.) The quantities of interest are then calculated on every ensemble configuration, and their mean values 
are obtained for each bin. Finally, the ensemble average is computed as the mean value of the bin averages, and its 
error is estimated as the standard deviation among the bin averages (if not noted otherwise). (On an Intel Core 2 
Quad Q6700 processor at 2.67 Ghz a thermal Markov step takes about ~ 7 seconds and a thermal sweep 2Ntm ^ 
1.89 hours.) 

IV. RESULTS AND DISCUSSION 

In the following subsections we analyze the physics content of dimeron ensembles which were numerically generated 
at the five squared coupling values — {1, 25, 100, 1000, oo} according to the procedure outlined in Sec. IIIII 

A. Dimeron dissociation as a function of the gauge coupling 

One of the essential features of the CDG mechanism is that with increasing coupling instantons are supposed to 
gradually dissociate into dimerons and to finally break up into their two meron partners. This process is suggested to 
be driven by a competition between the attraction among the regularized meron centers and their position entropy 
which increases with g^ (since the latter plays the role of temperature in the classical statistical analog ensemble). 
The competition is effective because both energy and entropy depend logarithmically on the distance 2\a\ between 
the meron centers. In fact, this distance will play a key role in characterizing both the structure of the ensembles' 
dimeron constituents and the interactions among them. The latter depend sensitively on the dimerons' color dipole 
moment which is of O (a^) . 

Before the dimerons break up completely, they should already have effectively released their two meron centers 
which then become the dynamically active degrees of freedom. In fact, this is expected to signal the onset of a phase 
transition in a finite system like ours where a complete |a| — > oo dissociation is prevented by the boundary. With 
their slowly decaying and thus strongly overlapping long-distance tails (^ l/^^) the essentially independent merons 
may then sufficiently disorder the vacuum to generate linear quark confinement. In the above sense, our dimeron 
configurations thus approach confining meron ensembles of the type studied in Refs. [3, E^- One should keep in 
mind, however, that our dimeron superposition ansatz (jl3p is not rich enough to describe ensembles of completely 
independent merons since it does not allow to untie the rigid color locking between the meron partners. 

Guided by the above considerations we are thus led to study the coupling dependence of the average dimeron 
dissociation 2(|a|) in our ensembles. These average inter-meron distances, which contribute to the typical dimeron 
size scale ~ 2 (p -I- (|a|)), are plotted as a function of the square coupling g^ in Fig. [5] and listed in Table HIl As 
expected, the average dimeron dissociation 2 (|a|) initially increases rather strongly with g^. For our two largest 
square couplings g^ = 1000 and oo, on the other hand, its value approaches saturation since it becomes comparable 
to the linear extent of the ensemble box. Hence 2 (|a|) will be underestimated in this coupling region. Moreover, large 
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dimerons close to the boundary are then increasingly forced to align themselves with the boundaries, and the largest 
dimerous accumulate along the box diagonals. Below we will nevertheless find indications for the meron centers (which 
are color zero-poles) to gradually replace the dimerons (i.e. color dipoles) as the dynamically most relevant degrees 
of freedom when the coupling increases, as alluded to above. 





2(|a|) 


1 


0.267 ± 0.003 


25 


0.567 ±0.001 


100 


0.893 ± 0.003 


1000 


2.239 ±0.014 


oo 


3.061 ±0.016 



TABLE II: The average inter- meron distance 2(|a|) for five values of the square coupling. 



B. Spacetime structure of the dimeron configurations 

We begin to explore the multi-dimeron physics of our ensembles by analyzing the spacetime structure of a typical 
member configuration. To this end we survey the action density of this field configuration and its amount of local 
(anti-) selfduality on different hyperplanes of the simulation volume. More specifically, we display these quantities 
as two-dimensional density plots in two section planes through the core volume which are chosen parallel to the 
a;4 — xi plane. One of them is separated from the — xi plane by the distances {x2, x^) = (0.4, 0.4) in the remaining 
directions, i.e. it lies rather close to the boundary in these transverse directions. The other one is positioned close to 
the center in the X3 direction, at {x2, X3) — (0.4, 1.2). (The origin of our coordinate system coincides with a corner of 
the core box. The spacial axes xi^2,3 range from to 2.5 and the X4 axis from to 4. Hence the box center is located 
at {xi,X2,X3, Xi) = (1.25, 1.25, 1.25, 2), cf.. Sec. MM) 

Since the detection and control of boundary effects is a recurrent issue when dealing with dimeron ensembles, we 
draw all plots for a 5^ = 1 configuration where the impact of the boundary is maximal (at least up to ~ 1000) and 
thus best analyzable. Such weak-coupling configurations with their small entropy are of additional interest because 
they are most likely to approximate semi-classical fields. In fact, at not too high densities (and not too large meron 
size regulators) such fields would be dominated by rather isolated and strongly contracted singular-gauge dimerons. 
These dilute superpositions of almost instanton-like solutions (cf. Sec. IIV A[) indeed approximate semiclassical systems 
quite similar to those studied in instanton vacuum models [23| . The following plots are designed to check how far our 

= 1 configuration resembles such semi-classical fields. 

Although our numerically generated dimeron field configurations contain redundant, i.e. gauge- dependent informa- 
tion without impact on observables, it is of technical interest to understand their spacetime structure since they lay 
the foundation for all our ensuing work. We have therefore examined typical dimeron ensemble configurations on the 
above set of hypersurfaces and found the fields to be fairly smooth. More importantly, all prominent spacial features 
of the gauge fields components were found to be closely mirrored in the gauge-invariant densities to be discussed 
below. In particular, we have found no evidence for the build-up of gauge-dependent peaks (potentially approximat- 
ing gauge singularities f45'|) which could adversely affect the simulation behavior even though they are invisible in 
gauge-invariant quantities. Hence we can refrain from plotting selected components of the gauge field itself. 

Instead, we turn to the gauge-invariant action density s (x) of our example configuration which is plotted in the 
two mentioned planes through the core volume in the two left panels of Fig. [6l The graphs indicate that the action 
density is indeed rather smooth, with the exception of a few dilute peaks. All these peaks show a slightly nonspherical 
shape and an extension of about 0.5. Now we recall from Table HIl that for g^ — 1 the average distance between the 
(anti-) meron centers of the (anti-) dimerons is 2 (|a|) ~ 0.27. This is about half of the regularized meron size p ~ 0.55 
and suggests to identify the peaks with the two strongly overlapping meron centers of the regularized dimerons. The 
average number of peaks is indeed consistent with the pseudoparticle density of the configuration. Moreover, the 
somewhat elongated action density of the peaks finds a natural explanation in the small but finite separation between 
the meron centers. Finally, in the more central {x2,X3) = (0.4,1.2) plane the peak density is larger than in the 
{x2,X3) = (0.4,0.4) plane which lies closer to the boundary. This may be a reflection of the strongly reduced average 
action density near the boundary (cf. Sec. IIIIDI and Fig. [3]). Statistical fluctuations are too large to substantiate this 
conjecture, however, as indicated by the fact that no such dilution is recognizable in Fig. [5] close to the boundaries in 
the xi,Xi directions. 

Another instructive property of the dimeron held configurations is their amount of local selfduality. This feature 
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characterizes the interplay between the Yang-Mills dynamics and the topological charge density q (x) (as defined in 
Eq. (|23|) ) and can be monitored by evaluating the expression 



which varies from —1 at positions where the field is selfdual to +1 where it is anti-selfdual. Since for = 1 the 
strongly overlapping, regularized meron centers render the (anti-) dimerous almost (anti-) instanton-like, one expects 
that the peaks are approximately (anti-) selfdual while the surrounding regions, dominated by overlapping tails, are 
neither. The plots of R {x) in the right panels of Fig. |6] confirm this expectation and thus allow to associate the 
action density peaks with either dimerous or anti-dimerons. Due to the one-to-one correspondence between the peaks 
in s {x) and R (x) the above comment on boundary effects applies here as well. 



We now proceed to a more quantitative analysis of the dimeron ensemble structure and its dependence on the 
square coupling . In the present section we look for patterns in the topological charge density which indicate short- 
and medium-range order. More specifically, we pick a pseudoparticle in a given ensemble configuration and measure 
the average density of the surrounding jDseudoparticles with either equal or opposite topological charge, as a function 
of their distance from the selected one 75] . We then repeat this procedure for all pseudoparticles in the configuration 
and average over the obtained density profiles, again for equal and opposite topological charge separately. After finally 
averaging over all configurations of the = 1 ensemble we end up with the two distance profiles pj^ plotted in the 
left panel of Fig. [71 The analogous procedure, but with the pseudoparticles replaced by individual meron centers, 
yields the profiles shown in the right panel. The integral of these densities over the full ensemble is normalized 
to one. (The smallest and largest distances are excluded in both figures since they correspond to tiny shell volumes 
(inside the ensemble volume) in which the densities cannot be reliably estimated.) Figures [8] and [9] contain the same 
profiles as Fig. [71 but for g"^ — 25 and g"^ = 100. 

The averaged radial density profiles reveal an intriguing amount of structure in the pseudoparticle distribution 
and in its meron substructure. First of all, the left panels of Figs. [71 ^ [H show a depletion of both dimeron and 
anti-dimeron densities in the overlap region with the reference (anti-)dimeron. Hence they provide direct evidence 
for a strong short-distance repulsion between dimerous of any topological charge. This repulsion is sensitive to the 
relative color arrangement between neighboring pseudoparticles (cf. Sec. IIVD|) and may at least partly be caused by 
our limited field configuration space. A similar repulsive core (with a logarithmic distance dependence) shows up in 
superpositions of instantons and anti-instantons in singular gauge (46j . On a practical level, it helps to avoid clustering 
among the pseudoparticles and promotes smoother and more semiclassical ensemble configurations. Evidence for the 
latter was already encountered in the density plots of Sec. IIVBI 

At larger distances, a remarkable medium-range order emerges among the pseudoparticles for g^ = 1 (left panel 
of Fig. [7]). In fact, about one meron size p ~ 0.55 (or the slightly larger average dimeron size p + (|a|)) from the 
fixed reference particle one finds an enhanced (depleted) density of pseudoparticles with opposite (equal) topological 
charge. At about 2p this structure is inverted and attenuated, i.e. the density of pseudoparticles with equal (opposite) 
charge is weakly enhanced (depleted). A further, even weaker inversion of the densities is discernible at distances 
~ 3/9, while from d > 3.5/9 both dimeron and antidimeron densities remain within errors equal to those of a random 
distribution. The almost periodic density oscillations over three consecutive layers indicate a pronounced mid-range 
order among the dimerous. In fact, the emerging shell structure resembles Debye-type screening clouds and indicates 
the existence of attractive short-distance correlations between dimerous and anti-dimerons (tg} . 

The above screening behavior should be enhanced at 5^ = 1 where the entropy is lowest and the field configurations 
therefore most strongly ordered. This can indeed be seen in our results. While for g"^ = 1 the third shells are clearly 
visible in Fig. [71 (although less pronounced than the first two), they essentially disappear for g^ > 25 (cf. Figs. [HI and 
[9]). Moreover, the dimeron densities in the left panels of Figs. [HI and [9l show a weaker first shell at somewhat larger 
distances (reflecting the growing average size of the dimerous), which now slightly favors pseudoparticles of equal 
topological charge. A hint of a second shell with inverted topological charge remains recognizable as well. Hence 
at stronger coupling and over typical nearest-neighbor distances the (anti-) dimerous show a tendency to surround 
themselves with (anti-) dimerous. This behavior could be another indication for the with growing dimeron dissociation 
increasing role of the meron centers as the dynamically relevant degrees of freedom. 

In order to test the latter interpretation, we show in the right panels of Figs. [71 -[HI the analogous density profiles 
Pj^i for individual meron centers (without regard for the dimerons which they are part of). At the smallest distances 




(18) 



C. Topological charge distribution 
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d <^ p, i.e. in the immediate overlap region with the fixed reference meron, the nearest-neighbor meron is most 
hkely its partner from the same dimeron, and hence has the same topological charge. This explains the enhanced 
(suppressed) density of equally (oppositely) charged merons for d <^ p. The enhancement is maximal for 5^ = 1 (cf. 
Fig. [7|) where the mean inter-meron distance 2 (|a|) is minimal, and decreases due the growing (|a|) with increasing 
. Outside of the immediate overlap region, on the other hand, i.e. for distances d > p, the meron density profiles 
in the right panels of Figs. [7] - 19] follow those of the corresponding dimeron densities (left panels) rather closely. This 
is expected because for d > p the encountered merons belong more likely to different dimerous. Since for — I 
with 2 (|a|) ~ 0.27 ~ p/2 the meron partners of the dimerous overlap almost completely, furthermore, their densities 
outside the immediate neighborhood of the reference meron should match those of the dimerous most closely, as is 
indeed the case. 

We have already pointed out that the strongest intermediate-range order among the dimerous exists in the g^ = 1 
ensemble with its particularly low entropy. For a more systematic analysis of the changes in this behavior with 
increasing coupling we now proceed to the investigation of global ensemble properties. (Their smaller statistical error 
simplifies the study of the g"^ dependence.) We first consider the fraction /^jj of dimerous in a given configuration 
whose nearest neighbor has opposite topological charge f77'|. The ensemble-averaged probability (/^jj) is plotted in 
Fig. [TU]for our five g^ values between one and infinity. (The vertical line indicates a break in the scale of the abscissa 
which allows to include the strong-coupling limit g^ = 00.) For g^ — 1 one reads off (/^q) — 87%. This large fraction 
confirms the strong preference of the dimerous to surround themselves with screening clouds consisting mostly of their 
anti-pseudoparticles. Already at g^ = 25 the value of (/qq) has diminished by half, however, to about 43%. This is a 
clear indication for the with increasing coupling growing disorder in the field configurations. It manifests itself not the 
least in the stronger dissociation of the dimerous (2 (|a|) ~ 0.57 ~ p at = 25, cf. Table HH which gradually replaces 
the interactions between dimeron centers by interactions between the increasingly independent meron centers. The 
fact that (/qd) remains within errors around 45% for g^ = 100 and 1000 signals the slight preference for equally 
charged pseudoparticle neighbors at stronger couplings, as already encountered in Figs. [8] and HI At these couplings 
the dimerous are so far dissociated, however, that correlations between their centers do probably no longer characterize 
the main interactions which they experience. In any case, the non-interacting value of 50% for (/qq) is attained only 
in the strong-coupling limit, i.e. in the random ensemble with g'^ = 00 (notwithstanding additional boundary effects 
which come into play for maximally dissociated dimerous, cf. Sec. IIII Pp . 

Further insight into the topological charge distribution of dimeron ensembles and its coupling dependence can be 
obtained from the average distances (d) between the pseudoparticle centers and their nearest neighbors with either 
equal or opposite topological charge. The nearest-dimeron distances are drawn in the left panel of Fig. [11] as a 
function of the square coupling g^. (They are evaluated in the same spherical shells as in Fig. [TOl) For g^ = 1 
one finds the average distance (d) ~ 0.71 between oppositely charged neighbors to be about 25% smaller than that 
between equally charged ones {{d) ~ 0.93). This is a reflection of the screening clouds found above. After increasing 
the square coupling to g^ = 25, on the other hand, the two distances have become almost identical. The nearest 
pseudoparticles of equal charge are now somewhat closer, as expected from the drop of (/djj) below 50%. (The 
average distance to the next dimeron of any charge remains relatively constant, though.) This seems to be another 
consequence of the larger inter-meron separation (2 \a\) ~ 0.6 — 3.1 for = 25 — 00 and the increasing dynamical 
importance of the meron centers. (For close-to-maximally dissociated dimerous the mentioned boundary effects will 
distort the nearest-neighbor distances as well, cf. Sec. IIII Dl ) 

Above we have found several pieces of evidence for the merons to gradually turn into the dynamically relevant 
degrees of freedom when the coupling increases. In order to explore this issue from yet another angle, we further 
computed the average distances between nearest-neighbor meron centers, as plotted in the right panel of Fig. [11] In 
the g^ — 1 ensemble with its highly contracted pseudoparticles, nearest-neighbor merons belong most likely to the 
same dimeron and thus carry identical topological charge. As a consequence, the average separation between equally 
charged neighbor merons is almost the same as the distance 2 (|a|) ~ 0.27 between partner merons. Oppositely charged 
neighbor merons, on the other hand, are more than twice as far separated, i.e. about as far as the average distance 
(d) ~ 0.59 between nearest-neighbor dimerous of opposite charge (cf. left panel of Fig. [TT|) . Again, this picture 
changes significantly with increasing coupling. As the dimerous dissociate farther, the average distance (d) between 
equally charged nearest-neighbor merons must increase as well. However, it does so more slowly than the average 
separation 2 (|a|) of the meron partners. This provides additional evidence for the interaction between individual 
merons belonging to different dimerous to increasingly determine the ensemble properties. At g^ = 25 the average 
distance between oppositely charged nearest-neighbor merons has not yet decreased much, on the other hand. The 
main drop happens between g^ = 25 and 100. From about g^ = 1000 the distances (d) between equally and oppositely 
charged merons conincide and slightly increase together on the approach to the strongly-coupled random ensemble. 
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D. Color correlations 



In the previous subsections we have found diverse spacial correlations in dimeron ensemble configurations and 
studied their coupling dependence. It remains to explore correlations between the color orientations of the dimerous 
to which we turn in the present section. There are several reasons for expecting these correlations to be significant. 
First, the "hedgehog" -type coupling between the spacetime and color dependence of the individual dimerous (cf. Eq. 
(O) suggests that spacial and color correlations will be at least partly linked in the ensemble configurations as well. 
In fact, we have encountered a prototype of such correlations when studying the color-orientation dependence of the 
classical interaction energy between two isolated dimerous. In weakly coupled ensembles, furthermore, the almost 
completely contracted dimerous will experience interactions and correlations similar to those found in singular-gauge 
instanton ensembles (see below). In the strong-coupling regime, on the other hand, where the dimerous are dissociated 
and single merons with their long-range tails become dominant, one may expect color correlations of similar strength 
and importance as in the meron and regular-gauge instanton ensembles of Refs. [l3. fioj. 

In order to characterize the SU(2) color orientation between two neighboring dimerous with color coordinates cp 
and Co (cf. Eq. ([H])), we introduce the angle 



a 2 arcsin 



(19) 



(with a e [Oi 71"]) under which half of the geodesic distance between the two corresponding points on the group manifold 
S'^ (of unit radius) appears from the center. We then find for every pseudoparticle in a given ensemble configuration 
the nearest neighbors with equal and opposite topological charges, separately calculate the angles (fT9|) between them 
and finally average over the configuration (excluding double-counting of pairs) and over the — 1 ensemble. 

The resulting (a) distributions are divided into 20 bins and plotted in Fig. [T2]for both equally and oppositely charged 
neighbors. For comparison, we also show the angle distribution of the non-interacting = oo random ensemble. All 
three distributions are relatively broadly peaked at (a) = 7r/2. This happens even in the random ensemble, which 
indicates that this peak position is statistically favored. Indeed, a — 7t/2 corresponds to the equator of and thus 
to the a value with the maximal number of color orientations. The main lesson of Fig. [T^] is, however, that equally 
and oppositely charged dimeron neighbors on average prefer remarkably different relative color orientations. While 
the distribution of the mean color angle (ajjjj) between nearest neighbors of opposite Q is essentially consistent with 
a random distribution, the peak value of the {ajjjj qq) distribution for neighbors of equal Q is dynamically enhanced 
by about a factor of two. 

These results can be rather directly understood by recalling that the interactions between the strongly contracted 
dimerous in the rather dilute g^ = 1 ensemble (with jd) > p, cf. Fig. [TT|) are similar to those among instantons in 
singular-gauge instanton-antiinstanton superpositions |23l . |46[ . Between instantons of equal topological charge these 
classical interactions ASjj jj are repulsive for any relative color orientation. More specifically, at distances d > p one 
has Ei 
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Since the average distance between DD and DD pairs in our g"^ = 1 ensemble is {djjjj Qo) ~ 0.93 ^ 2p (cf. Fig. 
[TT|) . the leading term in Eq. (^0]) should reasonably well approximate ASddi,^). The resulting a distribution is thus 
symmetric around 7r/2 and the individual repulsion is minimal (maximal) at a = 7r/2 (a = 0, tt). Hence the average 
repulsion in the ensemble is reduced by more strongly populating the a ~ 7r/2 color orientations. This explains the 
enhanced peak around (a^^ = 7r/2 in Fig. [T^l 

The instanton-anti-instanton interactions in singular-gauge instanton superpositions, on the other hand, contain 
the for d > p leading dipole-dipole interaction 
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where the unit vector d^ points from the instanton center to the anti-instanton center. For spacial and color orienta- 



tions with 
attractive. 



^fidf2 



IS 



< 1/2, which are predominant in the g^ — 1 equilibrium ensemble, the leading term in Eq. (|2T 

In fact, indirect evidence for this attraction was already deduced from our dimeron distributions in Sec. 
IIV CI All remaining // interactions are repulsive, on the other hand, as those between equally charged neighbors. 
In contrast to the latter case, however, the impact of AS^jj on the {ctui)) distribution is less transparent since the 
strength of the DD attraction is not determined by a. Moreover, the relatively small average dimeron-antidimeron 
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distance {dDjj) ~ 0.7 ^ 3p/2 for = 1 (cf. Fig. [TT|) indicates that the repulsive terms of O (p^/c?^) remain relevant 
as well. Therefore it seems likely that this repulsion, together with the overall short-distance repulsion (cf. Sec. IIVCP 
and the ensemble entropy, on average compensates at least part of the leading attraction. The resulting averaged 
DD correlations should thus be considerably weaker than those in the DD and DD channels. As a consequence, the 
{ajjjj) distribution will be close to that of a random ensemble, and this is indeed what we observe in Fig. 1121 

The above reasoning focused on two-pseudoparticle forces. This is legitimate because many-body interactions in 
the rather dilute = 1 ensemble are suppressed by the small packing fraction of the contracted dimerous and 
by their reduced overlap. Moreover, the emerging parallels with the color correlations in instanton-anti-instanton 
superpositions reinforce our premise that the properties of weakly coupled dimeron ensembles indeed approach those 
of instanton liquid models. As already alluded to, the above arguments further suggest that the smaller average 
distance (djjjj) < {dj^j^ Qo) between oppositely charged dimerous and the related screening-cloud arrangement in 
the = 1 ensemble (cf. Sec. IIV C|) are mainly caused by attractive color-dipole interactions. With growing coupling 
and entropy, however, the relative impact of the potentials ([20| and ([21]) on the free energy and on the ensemble 
structure decreases. This holds in particular for the DD attraction which has the longest range. Our finding of 
{^-dd) {d-DD dd) fo'" 9^ ^ 25 (cf. Fig. [TT|) indeed indicates that the common repulsion in both DD and DD,DD 
channels and the growing entropy begin to dominate at larger couplings. One would thus expect the peak in the 
{'^dd.dd) distribution to broaden with increasing coupling until the random distribution is reached for oo. 
The {olud) distribution, on the other hand, may become somewhat more sharpely peaked at larger g'^ when the 
compensating impact of the attraction subsides. 

In order to test these expectations and to she d further light on the coupling dependence of the color correlations, we 
have evaluated the standard deviation yj A {a) of the (a) distribution from its {g independent) mean value (a) = 7r/2. 
In Fig. [13] we plot A (a), again separately for equally and oppositely charged neighbor dimerous, at five g^ 
values. For — 1 the standard deviation of the {ctjjjj qq) distribution reaches only about half of that of the 
(ajji)) distribution, in agreement with its higher and narrower peak in Fig. 1121 Also expected from Fig. [12] is 
that the g^ = 1 value of -^Z A {ajjjj) remains close to the random ensemble value at = oo. In fact, the width 
of the distribution decreases only little with growing coupling, as anticipated above, until {ajjjj) becomes 

randomly distributed. The width of the {ctDD oo) distribution, on the other hand, grows rather strongly with g^ 
on its approach to the random value, again confirming our above expectations. Nevertheless, even at g^ = 1000 the 
{oiod,dd) distribution remains significantly more peaked than the random distribution. Hence some of the stronger 
average repulsion between nearest-neighbor dimerous of equal topological charge seems to remain influential even at 
these rather large couplings (see also Sec. lIVCp . 

Reflecting upon the above analysis one may wonder to what extent pertinent features of the a distributions are 
obscured by taking the configuration and ensemble averages. Indeed, the spacial averaging and the interplay between 
spacial and color correlations may wash out interesting local features of the relative color ordering. However, one 
would expect this type of leveling to be weakest in the g^ — \ ensemble where the almost spherical color distribution 
of the dimerous is most concentrated and the impact of their spacetime orientation consequently minimized. Other 
qualitative effects of the stronger intermediate-range order among neighboring dimerous aX g^ — 1 may therefore also 
be robust enough to survive the averaging procedure. 



E. Topological susceptibility 

The topological susceptibility Xt characterizes several fundamental properties of the Yang-Mills vacuum. It governs, 
for instance, the dependence of the free energy on the CP violating vacuum angle 9 around = and the mass of 
the rj' meson in quantum chromodynamics with a large number Nc of colors |5l| . Moreover, it contains information 
on the interplay between the Yang-Mills dynamics and the topology of the gauge group which was conjectured long 
ago to be the origin of linear quark confinement |2H . 

In a finite spacetime volume V , the topological susceptibility is defined as 

(which may contain contact terms [52l| in the p — and a — > limits) where Q {V) is the topological charge of the 
gauge field in i.e. 

Q(^)^ d^xg (x) with g (x) = -^tr {F^.i?,,.} (23) 

{F^i, = e^i,ai3Fafj/2). In the context of our vacuum description by topologically active constituents, x is of particular 
interest because it quantifies the strength of topological charge fluctuations around {Q) = 0. The latter is enforced 
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by CP symmetry which is sufHciently manifest in our ensembles because the configurations (jl3|) contain a very nearly 
equal number of dimerons and anti-dimerons j?^ . 

The topological susceptibility Yf = liniy^oo X (^) of SU(2) Yang-Mills theory is known from several independent 
high-precision lattice simulations [l^ll^l (with their scale set by prescribing the value of the physical string tension). At 
large Nc, furthermore, the value of Xt can be obtained either from lattice extrapolations [54-,56.] or from experimentally 
measured properties of the lightest pseudoscalar meson nonet. Indeed, to leading order in 1/Nc the Witten-Veneziano 
mechanism predicts [5l| 

Xt,su(iv.^oo) ^ ^ H, +ml~ 2ml) ^ (180 MeV)^ (24) 

(where /tt is the pion decay constant, Nf is the number of flavors and m^, m^' , rriK are the 77, rj' and K meson masses) 
and thereby relates Xt su(oo) to the part of the rj' mass which does not originate from the strange-quark mass. The 
above values of Xt '^i^ serve as benchmarks for comparison with our results and for scale-setting purposes. 

To evaluate x in our dimeron ensembles, we first compute the topological charge Q {V) of each gauge-field con- 
figuration. More specifically, we integrate the topological charge density, obtained analytically from Eqs. (IT5|) and 
(1^51) . on a special grid [t^ whose extent is varied in small steps to tune through the desired range of volumina V. 
The correlation between discretization errors in the two Q factors of x (which arise from neglecting fluctuations with 
wavelengths below the grid-point distance in the coarse-grained density q) is reduced by evaluating the second factor 
on a modified grid. The latter is obtained by replacing each of the original grid-point positions with a randomly 
chosen one inside a surrounding volume determined by the inverse grid-point density. Finally, the ensemble average 
is carried out according to Eq. p2|). 

Translational invariance implies that x (^) becomes volume-independent in the thermodynamic limit and that 
Xt = limv-^oo X (^) exists (after appropriate renormalization). In fact, this V independence will develop already in 
finite volumes whose linear dimensions ^ V^^^ sufficiently much exceed the correlation length of the topological charge 
density. In order to test how far this applies to our simulation volumes and permits extrapolations of our x (^) to 
V — >■ 00, we have calculated x (^) in a range of volumes in which boundary effects remain controllable. (Numerical 
efficiency is improved by varying only the temporal extent of the evaluation volumes. All three spacial dimensions of 
the boxes are kept fixed at 1.2 (chosen to limit boundary artefacts, cf. Sec. IIII Pp and remain centrally embedded.) 

The resulting volume dependence of the topological susceptibility is shown in Fig. [14] for = 1 and 100. In small 
volumes topological charge fiuctuations are suppressed and x (V) accordingly starts out close to zero. Towards our 
largest reliably accessible volumes, on the other hand, boundary effects may become relevant and begin to reduce 
X- (Hence Fig. [TJ] contains independent information on the extent to which the boundary breaks translational 
invariance.) Nevertheless, for both coupling values x (^) shows a rather broad maximum towards the end of the 
trustworthy V range. It is tempting to interpret these maxima as the onset of the expected saturation plateaus, and 
the corresponding x values as reliable approximations to the infinite- volume predictions Xt- However, from our data 
sets one cannot decide with confidence whether the maxima do not instead underestimate Xt by interpolating between 
the rise of x (V) at small V and its boundary-induced decline at large V. We shall therefore interpret the maximal 
values X (in this and the following sections a hat will indicate a dimensionful quantity given in our numerical units) 
more conservatively as lower bounds on the topological susceptibility. 
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a-x [fm] 




1 


0.75 ± 0.032 


0.076 ± 0.0032 


3. 864 ± 0.0793 


25 


0.83 ± 0.043 


0.084 ± 0.0045 


9.832 ± 0.7271 


100 


0.94 ± 0.032 


0.095 ± 0.0033 


9.340 ± 0.4518 


1000 


2.42 ± 0.083 


0.246 ± 0.0096 


3.363 ±0.1051 


00 


10.35 ±0.503 


1.043 ± 0.0508 


0.1607 ± 0.0008514 



TABLE III: (Lower bounds on) the fourth root of the topological susceptibility x in numerical units for five values of the 
squared gauge coupling. The last columns contain the grid constant and the string tension a in physical units as obtained 
when setting the scale with the SU(2) Yang-Mills value x^''* = 0.99 fm"\ 

In Table IIIII the corresponding values of x'^'^ are listed with statistical (jackknife) error estimates for our five 
standard values. On a cautionary note, we recall that the g^ = 1 dimeron ensemble thermalizes very slowly and 
that the autocorrelation for the topological charge becomes correspondingly large. The error of x^^^ as quoted in 
Table IIIII does not take this strong autocorrelation into account and is therefore probably underestimated at least for 
g^ — 1. With this caveat in mind, we observe that x'^'^ stays inside our error estimates for 25 < g^ < 100 practically 
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constant (which becomes even more obvious when expressing y^^'^ in physical units, cf. Table IV)) . This weak coupling 
dependence may be related to the similarly weak temperature dependence of the topological susceptibility in the 
confined phase of Yang-Mills theory up to the phase transition [s^. For larger our values grow too strongly, 
on the other hand, which is at least in part due to boundary artefacts. Indeed, from > 1000 the dimerous are on 
average as far dissociated as the boundary allows. As discussed in Sec. IIIIDl they then tend to accumulate along the 
box diagonals, with their centers concentrating in the middle of the box. These distorted field configurations break 
translational invariance almost everywhere in the simulation volume. 

Our numerical results for the topological susceptibility provide a first opportunity to set the distance scale a in 
our box as a function of g^. To this end we write = a (x/x)^^* in physical units, where a — 0.1 and x ^re given 
in our numerical units while x — (0.99)'* fm"** is the SU(2) Yang-Mills value obtained from the SU(2) lattice result 
X^/'*/(T^/^ — 0.483 ± 0.006 j54j in combination with the physical string tension a ~ 4.2 fm~^ [l[. The resulting values 
of in fm are listed in Table Hill With their help one may convert our data e.g. for the string tension a (cf. Table 
IIVP via a = {a/a-^)^ u into physical units, as done in the last column of Table lllll 

A quantitative comparison of our results with those of other approaches in physical units will be postponed to Sec. 
IIV Gl For a first orientation, however, one may compare our "raw data" for x with their counterparts from the regular- 
gauge instanton and meron ensembles of Ref. ^] (which adopted = 32 and a size regulator p < 0.1) even though 
the physical distance scales might not be straightforwardly related. (In the ensembles of Ref. JJJ about 90% of the 
topological susceptibility is generated by the peaks in the topological charge density associated with individual merons 
or regular-gauge instantons, incidentally; the remaining, smoother background field contributes just the remaining 
10%.) According to Table Hill our ensembles with g^ — 2h — 100 give x^/'^ ~ 0.8 — 0.9 which is of the same order 
as the meron-ensemble result x^^^ = 0.77 14 1 for a comparable meron density in physical units, corresponding to 
-^M+M = 100 (sse below). In fact, at the above couplings g^ = 25 and 100 our dimeron density reaches about half of 
Ref. [14|]'s meron density and therefore yields a similar density of meron constituents. (The regular-gauge instanton 
ensemble with an instanton density roughly equal to our dimeron density (corresponding to Nj^j = 100), on the other 

hand, generates the significantly larger value x^^^ — 1-37 14].) Within errors the topological susceptibilities of meron 
and dimeron ensembles in numerical units are therefore consistent. As a consequence, setting scales by imposing a 
common "physical" value for x will lead to comparable physical distance scales. 



F. Wilson loops, static-quark potentials and string tension 

The search for the conjectured transition of dimeron ensembles into confining meron ensembles is one of the central 
objectives of our study (cf. Sec. Our first quantitative evidence for the actual development of a meron-dominated 
phase was discussed in Sec. IIVBI where we found the dimerous to dissociate into their progressively independent 
meron constituents when the coupling increases. We are now proceeding to the crucial question whether the gradually 
liberated merons indeed generate an area law for large Wilson loops and thus linear confinement with a finite string 
tension. More specifically, we are going to evaluate ensemble averages (i.e. vacuum expectation values) of path-ordered 
Wilson loops 

W [C {R, T)] =tr{Vexpi(l> A^dx^, | (25) 

along closed, rectangular paths C. The latter correspond to the worldlines of fundamentally colored, static quark and 
antiquark sources, separated by a spacial distance R and evolving over the Euclidean time T > R. The static energy 
required to introduce and separate these sources in the vacuum is then given by the potential 

V (R) = - \im ^In {W[CiR,T)]). (26) 

Hence a linear rise of V (i?) for sufficiently large R signals linear quark confinement by color-flux tube formation. 

Since we are particularly interested in the behavior of V (R) at the largest reliably calculable R values, our strategy 
for computing W is designed to minimize the impact of boundary effects (as identified in Sec. IIII Pp while still 
making efficient use of the ensemble configurations. This is achieved by restricting the evaluation volume to several 
subvolumina centered in the core box of Sec. IIII Al Artefacts of the boundary are reduced by limiting the extent of 
these volumina to 2.0 in the time direction, to 1.0 in the spacial R direction and to 0.5 in the two remaining spacial 
directions. The volumina generated by permutations of the spacial directions are included as well. As a consequence, 
fields in regions in which we diagnosed the main violations of translational invariance are excluded from the evaluation. 
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Next, we compute the oriented links U [I {x —S' x)] = trV expi Jk^^^^,) Apdxp which connect neighboring grid points 



and in the evaluation volume by straight lines I {x 



All rectangular Wilson loops with aspect ratio 



R/T =1/2 are then constructed from combinations of these color parallel transporters. (For small enough loops we 
exploit Euclidean symmetry to assign the time direction to the longer side, irrespectively of its orientation in the 
core box.) In order to obtain quantitative information about the impact of boundary and resolution errors on the 
reliability of the link calculation, we perform the latter both (i) on the full subvolumina described above, with the R 
values progressing in steps of length 0.1, and (ii) just on the central hyperplanes of these volumina, spanned by the 
longer spacial direction (with length 1.0) and the temporal direction, in which the R values progress with step size 
0.025. 

The main challenge for our finite grid resolution arises from gauge fields along a link which develop unusually large 
gradients by approaching a meron center, as already alluded to in Sec. IIIIDl Since this turns out to happen relatively 
rarely, we decided to dynamically adapt the resolution of the grid on which the gauge fields are sampled along the 
link. The refined grid then increases the accuracy of the numerical integration in the link exponent. We start with 
a maximal grid-point distance of 0.01, i.e. on a grid which is 10 times finer than the one on which the gauge-field 
configurations were generated. This distance is reduced by half if the Frobenius norm of the difference between a 
given link and the path-ordered product of two links over half the distance exceeds one per mille, i.e. if 



\U[l {x, 



X.J 



i)] -U[l {xp ^ 2:^+1/2)] U [I (x^+1/2 ^ Xp+i)] II > 0.001. 



(27) 



The criterion ((27)) is chosen such that about one third of the links in an average Wilson loop has to be updated at 

least once. (Two refinement updates of the subgrid for the same link are almost never encountered with our regulator 
= 0.3.) Finally, we calculate the average over all Wilson loops [s^ with the same edge length R and evaluate the 

logarithm In {W [C {R, 2R)]) of their ensemble average. 

In Fig. [15] the resulting data are plotted as a function of the loop area A = 2R^ for our five square couplings 
— 1, 25, 100, 1000 and 00. The figure includes the points obtained from both the full subvolumina and their central 

hyperplanes, which agree within errors 8l|. It further shows the best fits of the In (W) data to the function 



f{A)=uj + TP (A) - a A 



(28) 



where P (A) = 3V2A is the perimeter of the rectangular loops. The ansatz (f28l) may e.g. be motivated by a string 
model which takes small, UV- regularized quantum fluctuations into account [l4l53l- All fit results are collected in 
Table HVl together with their values per degree of freedom and the fit ranges A G [Amin,^max] in which the data 
are probably reliable. The values indicate that the fits indeed allow for a reliable determination of the parameters 
w, f and in particular of the string tension a. The central lesson of these results is that dimeron ensembles indeed 
generate finite string tensions which monotonically increase with g^. Our findings in Sees. IIV Al - [IV Dl explain this 
increase as a consequence of the dimerous' gradual release of their disordering meron tails. In addition, our results 
provide evidence for dimeron dissociation to trigger the transition from a non-confining regime to a confining phase. 
(In the thermodynamic limit the analogy with Kosterlitz-Thouless-type transitions suggests that the string tension 
will vanish exactly below some critical value g^. In finite volumes one expects a more gradual transition, however, as 
seen in our data.) 





Co 


f 


a 






-^max 


1 


-0.22 ±0.09 


0.27 ±0.07 


2.23 ±0.22 


0.22 


0.20 


1.60 


25 


-0.51 ±0.04 


0.81 ±0.04 


6.95 ±0.21 


0.32 


0.10 


0.90 


100 


-0.20 ±0.01 


0.55 ±0.02 


8.47 ±0.18 


0.44 


0.03 


0.51 


1000 


-0.02 ±0.04 


0.47 ±0.08 


20.08 ±0.66 


0.21 


0.03 


0.35 


00 


-0.63 ±0.09 


0.23 ±0.20 


17.57 ± 1.80 


1.06 


0.01 


0.35 



TABLE IV: Results for the parameters u, f and the string tension 0, obtained from the fit of f{A) to \n(yV){A) in the interval 
from Auiiu to ylmax. (The reduced chi-square value characterizes the statistical goodness of the fit and should not be confused 
with the square of the topological susceptibility.) 



We have additionally applied a somewhat complementary method of computing the string tension. In this approach 
the static quark-antiquark potential V (R) is extracted directly from the Wilson-loop behavior according to Eq. . 
This allows us to include loops of multiple aspect ratios, provides a useful cross check on our results for a in Table 
IIVI and makes it possible to obtain physical information from the small- i? behavior of V as well. To implement this 
approach we evaluate rectangular loops of all aspect ratios which fit into the subvolumina defined above. The value 
of F at a given R is then extracted from a fit of In (W) /T to a constant inside the T > R ranges in which In (W) /T 
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becomes approximately T independent. The resulting potentials are plotted in Fig. [16] for = 1,25 and 100. All 

V (R) curves indeed show an essentially linear rise for R > 0.4. (Due to growing boundary effects, the potentials 
cannot be reliably extracted beyond a coupling-dependent i?max ^ 0.7 — 1.) It is reassuring, furthermore, that fits 
to the linear region yield values for the string tension (tr (g^ = l) ~ 1.6, a {g^ ~ 25) ~ 6.3 and a [g^ — lOO) ~ 9.3) 
which are within 10 ~ 20% of those with better statistics given in Table HVl 

Before setting scales and comparing to results of other approaches in Sec. IIV G\ we may get a first idea of the 
quantitative significance of our Cj,t and especially a values by comparing them to their counterparts in the meron 
and regular-gauge instanton ensembles of Ref. [lj| (which adopted g^ — 32, p < 0.1 and Njj^j j^^^j^j = 500). The 
meron ensemble yields oj = —0.5, f = 0.94 and a ~ 12.8 which are overall closest to our — 25 and 100 results. 
However, the string tension in the meron ensemble exceeds ours in the g^ — 100 dimeron ensemble by about 50%. 
This is probably because our dimerous' constituent merons are not yet fully Hberated at ^ 25 - 100 (cf. Table 
|TT| and because their density is somewhat lower than the meron density in the Nj^_^_j^f = 500 ensemble of Ref. ■ 
In any case, the "raw" data suggest that the confinement properties of g^ = 32 meron and g^ = 25 — 100 dimeron 
ensembles are at least qualitatively compatible. (Another indication may be that the string tensions in both meron 
and dimeron ensembles roughly double in the strongly-coupling limit, boundary effects notwithstanding. The results 
cj = —0.1, T — 0.58, a = 20.5 for the regular-gauge instanton ensemble 14], incidentally, turn out to reproduce ours 
from the g^ ~ 1000 dimeron ensemble up to a few percent.) 

The above considerations indicate that fixing the string tension at a common physical value will result in comparable 
confinement scales of both meron and dimeron ensembles in the g^ — 25 — 100 coupling region. This is reassuring 
because the linear rise of our V (R) should then rather directly match on to the linear potential of the meron ensemble 
even at i? > 1, where the confining flux tube is expected to develop fully but where we cannot reliably extract it from 
our present data set. This expectation is supported by fits of the meron-ensemble result for V (R) to a regularized 
string model in an extended R range including R > 1 [1J|. In fact, within errors those fits reproduce the string tension 
obtained from fitting In (W) to Eq. (pS)) in an A range with a smaller i?niax — \/^/2 than ours. This indicates that 
our string tension can indeed be reliably determined by fitting Eq. (|28p in a region of relatively small areas. It 
also adds considerably to the evidence for dimeron ensembles to provide a viable pathway towards confining meron 
ensembles. 

We close this section with a few comments on the small- i? behavior of V (i?) which is of physical interest in its 
own right. For R < 0.4 our potentials exhibit to a good approximation a quadratic R dependence which extends 
farthest for g^ = 1 and is little affected by boundary artefacts. In order to understand this behavior one should recall 
that such a quadratic small- i? potential ^3] is generated by sufficiently dilute instanton-antiinstanton ensembles (23 - 
[26j . This indicates that for i? < p our V {R) are dominated by contributions from little dissociated dimerous which 
behave essentially as singular-gauge instantons. (Table [ll] and the dimeron density estimates in Sec. IIV Gl suggest 
that our ensembles keep a substantial fraction of such instanton-like dimerous even at larger couplings.) For a more 
quantitative consistency check we recall that a dilute SU(2) instanton ensemble generates the short-distance behavior 

V {R) ~ 0.5p~'^i?^ [l^l, i.e. V {R) ~ 2.8 i?^ for our size parameter p = 0.55. From a quadratic fit to our g"^ ~ 1 
potential at distances R < 0.4, on the other hand, we obtain V (R) ~ 3.8 R^ which is indeed of the same order. 



G. Dimensionless amplitude ratios and results in physical units 

A quantitative comparison of the dimensionful amplitudes and observables calculated above with those of other 
approaches requires either to form dimensionless combinations of our results or to rewrite them in physical units. 
Both will be done in the present section. We start with the discussion of three dimensionless ratios which can be 
formed from our results and thus compared to other data without scale-setting ambiguities. 

Our first example is the ratio (s) /ct^ — (s) /tr^ of the ensemble-averaged Yang-Mills action density (s) (the hat 
again indicates numerical units) and the appropriate power of the string tension given in Table IIVI The density g^ (s) 
is (modulo renormalization issues) proportional to the gluon condensate (-f^) and therefore of considerable interest 
in its own right. Our results for g^ (s) can be obtained from the values computed in Sec. IIIIDI on diagonals through 
the simulation volume by taking the average values of the plateau fits in Figs. OandlH i.e. g^ (s) (g^ = l) ~ 210 and 
^2 (^2 ^ 25) ~ 425. With the resuhs for a in Table |IV] this yields (s) /ct^ (^2 i) ^ 94 and (s) /ct^ (g^ ^ 25) ~ 3. 
The latter value corresponds to the coupling range 5^ ~ 25 — 100 in which the dimeron ensembles were estimated 
above to best approximate the confining phase of the Yang-Mills vacuum. Hence this value should be compared to the 
SU(2) lattice value (s) /ct^ ~ 4.5 60], the QCD sum- rule values (s) / ~ 4 — 10 [11] (which correspond to A^c = 3 and 
contain quark corrections, however) and the value (s) /cr^ ~ 8 [83[ in the meron ensemble jl3|- (The regular-gauge 
instanton ensemble yields (s) /cr^ ~ 10 [l3|-)- Of course, all these values should be taken with a grain of salt since 
they involve large subtractions whose systematic error cannot be reliably estimated. In any case, within expected 
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errors our result for = 25 is compatible with the SU(2) lattice value and several QCD sum-rule estimates. 

A second useful dimensionless ratio to be assembled from our results is xl (ff^ (s)) • In fact, there has been a proposal 
for an approximate low-energy relation in Yang-Mills theory, x ~ (-F^) / (BGtt^TVc) [61|, which would fix this ratio as 
xl {g^ (s)) ^1/ (SStt^) ~ 3.07X 10~3 (for TV^ = 2 and with (s) = (F^) /4). A dilute-instanton-gas estimate similarly 
yields x/ (ff^ (s)) - 1/ (327r2) [§3. With the SU(2) Yang-Mills lattice rcsuh x^^'^/cr^^'^ = 0.483 ± 0.006 [H^ and the 
ratio quoted above one arrives at the similar value x/ id^ (*)) — 3.21 x 10^^ while the range of condensate results 
{F"^^ ~ 0.1 — 0.3 GeV' 49] from QCD sum-rule analyses (corresponding to A^'c = 3 and containing quark admixtures) 
together with the SU(2) Yang-MiUs value x^'^ - 195 MeV (see below) yields xl {g^ (s)) - (1-2 - 3.6) x 10"^. From 
the data for x in Table Hill and for (s) as given above, finally, our dimeron ensemble results are xl (ff^ (■s)) — 1-9 ^ 10~^ 
for = 1 and x/ {d^ (s)) ~ 1. 2 x 10~^ for g^ = 25, within the range obtained from the sum-rule values for the gluon 
condensate but smaller than the lattice value. 

The third dimensionless quantity which can be formed from our results is x^/^/u^/^. In contrast to the ratios 
discussed above, it has the additional benefit of not involving the rather unreliably known and renormalization-scale 
dependent expectation value of the action density. Our x^^"*/""^/^ values for the five different g'^ are listed in Table fVl 
The inverse ratio cr^^^/x^^'' is plotted for the four finite g'^ values in Fig. [171 For g^ between 1 and 1000 our values 
for x}^^ jo^l"^ lie in the range from 0.3 to 0.55. Hence they are indeed rather compatible with the SU(2) lattice value 
X^/'*/cr^/^ — 0.483 ± 0.006j5J]. Moreover, for g^ — 2h— 100 our results are practically identical to the meron-ensemble 
value x}l^ ja^l"^ ~ 0.31 [l^ at g^ = 32 (and somewhat smaller than the range 0.42 < x}^^ 1^^^"^ ^ 0.48 obtained from 
regular-gauge instanton ensembles [l3|). 

Since correlations among individual observables may be hidden in the dimensionless ratios discussed above, it is also 
important to directly compare dimensionful results in physical units with those of other approaches. For this purpose, 
we have considered two alternative scale-setting procedures. The first was introduced in Sec. IIV El and prescribes a 
standard value for the topological susceptibility in SU(2) Yang-Mills theory. This allows to rewrite our results for the 
string tension in physical units, as given in the last column of Table Hill For = 25 — 100 they are roughly twice as 
large as those obtained from Regge phenomenology (see below), while for g^ — 1000 our value is about 25% smaller. 
Continuing this trend, the string tension becomes unnaturally small in the strong-coupling limit, as a consequence of 
X growing too large due to boundary artefacts (cf. Sec. IIV El and below). 

In the following, we therefore adopt a more standard approach to scale setting which fixes the string tension at 
its physical value a ~ 4.2 fm^^ (estimated from the almost universal experimental Regge slopes l]). This amounts 

to writing our discretization scale in physical units as a = a (o•/^7)^/^ with a = 0.1 and a (cf. Table HV]) given m 
numerical units. The resulting values in Table [Vl show that a increases with g^ . Hence the resolution of short-distance 
features in the vacuum field population decreases moderately with the coupling. This tendency is enhanced by the 
gradual dissociation of the well-localized, singular-gauge dimerous into broader merons (cf. Sec. IIV A[) which contain 
less short-wavelength Fourier modes. As argued above, we expect our ensembles to best reproduce observables in the 
confining Yang-Mills vacuum phase for = 25 — 100. In contrast, our most weakly coupled ensemble with g^ = \ 
underestimates the string tension while for g^ ^ 100 the results increasingly suffer from boundary artefacts. In those 
ensembles, fixing the string tension at its physical value will therefore not yield dimensionful results of the expected 
magnitude. With this caveat in mind, we can now use a to express other dimensionful quantities in physical units (cf. 
Table El). 





1/4 / 1/2 


a [fm] 


pM 


n [fm^^] 


[MeV] 


1 


0.52 ±0.057 


0.073 ±0.004 


0.040 ± 0.002 


15.69 ±3.44 


209.02 ± 23.04 


25 


0.29 ±0.022 


0.129 ±0.002 


0.071 ±0.001 


1.71 ±0.099 


117.61 ±8.90 


100 


0.32 ±0.012 


0.142 ±0.002 


0.078 ± 0.001 


1.10 ±0.062 


127.24 ± 5.034 


1000 


0.54 ±0.033 


0.219 ±0.004 


0.121 ±0.002 


0.19 ±0.014 


218.66 ± 13.44 


oo 


2.46 ±0.13 


0.205 ±0.01 


0.113 ±0.006 


0.25 ± 0.049 


991.28 ±51.42 



TABLE V: Results in physical units set by the string tension a — 4.2fm 



An important dimensionful quantity is the meron-center size regulator p which furnishes the characteristic input 
scale of the dimeron ensembles (cf. Sec. Ill C[) . For 5^ = 25 — 100 and p — 0.55 we find p ~ 0.07 — 0.08 fm, i.e. about 
a quarter of the typical instanton size pj ^ 0.3 fm in (singular-gauge) instanton-liquid models 23|. This indicates 
that our regularization procedure deforms the classical dimeron solutions only moderately and should not impede a 
potentially semiclassical behavior of their superpositions. The average dissociation 2 (jaj) ^ (1 — 1.6) p of the dimerons 
at g^ — 25 — 100 increases the effective dimeron size, on the other hand, and brings it closer to typical instanton 
sizes. Finally, it is useful to note that /i ~ p~^ ~ 2.5 — 2.8 GeV approximately sets the renormalization scale of our 
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ensembles. (Recall that p acts as an UV regulator e.g. for the action density which diverges logarithmically when 
p — >■ 0.) These /i values are of the order of typical lattice scales, which may help to explain why our results for the 
scale-dependent action density are compatible with lattice values. 

We now turn to the dimeron density of our ensembles. Recalling the ^rid-point distribution from Sees. IIII Al and 
IIIIEl the four-volume of our ensemble box is given by 14ns = (1.15 x 25) x (1.15 x 40) a^. The corresponding values 
of the dimeron density n = [No + Njy) /Vcns — 487/T4ns in physical units are listed in Table |Vl In the preferred 
coupling region = 25 — 100 they are about 20-80% larger than the typical instanton densities nJ_^_J ~ 1 fm""* of 
instanton liquid models [23]. The physical significance of this result will be discussed below. Moreover, it is interesting 
to note that the density nj^j^fj ~ 3.2 fm""' of merons in the Nm — 200 ensemble of Ref. 141] is similar to the density 
of our meron constituents, i.e. approximately twice as large as our dimeron densities n ~ 1.2 — 1.8. This adds to our 
earlier evidence for an increasingly dynamical role of the largely independent meron constituents from the farthest 
dissociated dimerous. 

It remains to put our results for the topological susceptibility in the last column of Table |V] into a quantitative 
perspective. The dimeron ensemble predictions for x^^^ in the coupling region = 25 ~ 100 are about 50% smaller 
than the "physical" value x^^"^ = 195.03 MeV (i.e. the SU(2) lattice value with the scale set by cr = 4.2 fm~^, cf. 
Sec. HVE)) . Similarly smah values x^^'^ ~ 118 - 132 MeV were encountered in the meron ensembles of Ref. [T^ 
while regular-gauge instanton ensembles yield x^^^ — 162 — 190 MeV. Hence our relatively small x^^^ values are 
consistent with previous indications for the meron centers, which carry only half of the dimerous' topological charge, 
to increasingly determine the topological-charge fiuctuations for g^ > 25. (One should also keep in mind, however, 
that we consider our values for x as lower bounds, cf. Sec. IIVEI ) Another suggestive pattern in our data is the weak 
bare-coupling dependence of the x^/"* values for = 25 — 100. As already alluded to, this behavior ma.y reflect the 
well-known temperature independence of the topological susceptibility in the confined Yang-Mills phase |56j which in 
turn indicates that x^^^ is relatively insensitive to temperature-induced variations of the coupling gym (-^) '^P the 
critical temperature Tc- (For T > Tc, on the other hand, x drops very rapidly.) 

The unusually large value of x^^^ ii^ our g^ — oo ensemble was already noted above in the raw data (cf. Table 
imp . Of course this result raises suspicion. In a random distribution of merons, into which the dimerous dissociate 
for 5^ — > oo as far as our limited field configuration space and the boundary allows, one would expect the topological 
susceptibility to be an incoherent sum of single-pseudoparticle contributions and thus to have a value of the order of 
the pseudoparticle density. This expectation was verified explicitly in random singular-gauge instanton [23j . regular- 
gauge instanton and single-meron |14| ensembles j85j]. The fact that our for g^ oo strongly increasing x values 
fail to scale with the density is therefore probably another manifestation of the residual correlations between the 
dimerous' meron partners and the related — )■ oo boundary artefacts. We recall that in Sees. IIV Al and IIV El these 
artefacts were found to impede the development of a genuine random meron ensemble and to strongly contaminate 
the topology distributions. Fortunately, in the physically relevant coupling region g^ = 25 — 100 these effects are much 
weaker and under far better control. 

Our above results suggest a remarkably comprehensive role for dimerons in the Yang-Mills vacuum. A key feature of 
the emerging picture is the rather distinct division of labour between almost fully contracted, instanton-like dimerons 
and their far dissociated, meron-like counterparts. In dynamical equilibrium at intermediate, i.e. physical couplings 
both of these components coexist at fairly high densities [s^. In particular, our overall dimeron densities significantly 
exceed the instanton density in instanton liquid models. Hence the instanton-like dimeron component should be 
sufficiently populated to perform essentially all established tasks of instantons in the Yang-Mills vacuum ^23;, 231 ■ The 
strongly quadratic rise of our heavy-quark potentials at small interquark separations provides additional evidence in 
this direction (cf. Sec. lIVFp . The meron-like component, on the other hand, seems to be mainly responsible for the 
longest-distance vacuum physics including, most importantly, linear quark confinement as established in Sec. IIVFI 
and in single-meron ensembles. 

On a structural level, the above reinterpretation of the instanton component in the Yang-Mills vacuum as mainly 
consisting of hardly dissociated dimerons is compatible with the existing instanton phenomenology as well. Indeed, 
due to the destructive interference between their meron-center tails the instanton-like dimerons experience much 
weaker long-range interactions than the almost liberated merons of the meron-like component (cf. Sees. lIVBI - ITVDp . 
Similar to singular-gauge instantons in sufficiently dilute superpositions they overlap less with each other, furthermore, 
and thus retain more of the classical solutions' shapes and semiclassical behavior. Hence the instanton-like dimeron 
component indeed resembles a rather weakly interacting "liquid" . On the other hand, this raises the question why 
cooling studies, designed to detect semiclassical solutions in equilibrated lattice configurations, have reported clear 
evidence for instantons or calorons with non-trivial holonomy but none for dimerons [l6ll29j. Closer inspection reveals, 
however, that these findings are not in conflict with our above scenario. In fact, the finite-resolution algorithms used 
to identify these solutions by their shapes may have counted small dimerons as instantons. Moreover, under cooling 
the logarithmic attraction between the dimerous' (smoothed) meron centers can no longer be counterbalanced by 
their decreasing entropy. Hence the dimerons will shrink and at least partially coalesce into instantons [g^l as which 
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they are identified afterwards. 

As established in Sec. lIVFi the effectively released meron centers of the meron-like dimeron component develop 
the strong long-range correlations required to generate confinement. Although these color correlations produce a 
more complex long-distance structure than encountered in fully disordered random ensembles, they still generate 
enough entropy for a string tension of the physically expected magnitude to develop (cf. Sec. lIVFp . Moreover, it 
seems likely that our above scenario can be extended to an at least qualitatively realistic description of the QCD 
vacuum. The instanton-like component of dimeron ensembles, in particular, should adapt to the presence of light 
quarks in essentially the same manner as in instanton liquid models. This is because most instanton effects in the light- 
quark sector, including the induced quark interactions which spontaneously break chiral symmetry, originate from the 
characteristic quark zero-modes arising in the instanton background. Since these zero modes are of topological origin, 
their existence and bulk properties [6J| are little affected by moderate deformations of the instantons into dimerous. 
Hence the ensuing phenomenology should be similarly robust under such deformations. 



V. SUMMARY AND CONCLUSIONS 



We have studied the extent to which the Yang-Mills vacuum can be described by restricting its gauge-field content to 
superpositions of regularized SU(2) (anti-) dimeron solutions. Our set-up and analysis of the corresponding effective 
theory focused mainly on structural changes in dimeron ensembles as a function of the bare gauge coupling and 
traced their impact on topological and confinement properties. The underlying quantum superpositions of dimerous, 
with their dynamics governed by the Yang-Mills action, were generated without recourse to weak-coupling or low- 
density approximations by Monte-Carlo simulations in a wide range of bare coupling values. The localization of 
the individual dimerons was enhanced by transforming them into a singular gauge, which simultaneously improves 
numerical efficiency, the vacuum description at shorter distances and the compatibility with semiclassically motivated 
instanton liquid models. 

Our initial survey of dimeron ensemble properties and their coupling dependence concentrated on the spacial 
distribution of the dimerons, their local degree of selfduality and their meron substructure. To understand relevant 
aspects of the underlying inter-dimeron dynamics, we further studied short- and intermediate-range correlations 
between positions, topological charges and color orientations of neighboring dimerons and of their meron centers. 
The resulting structural patterns turned out to be multifaceted and revealing. A strong repulsion up to distances of 
the meron center sizes governs the short-range interactions between dimerons of both topological charges and at all 
couplings. In our minimally coupled and thus most ordered ensemble, we additionally find a striking (about sevenfold) 
preference for nearest-neighbor dimerons to carry opposite topological charges, indicating a strong intermediate-range 
attraction (repulsion) between dimerons of opposite (equal) topological charge. A less pronounced but still significant 
majority of next-to-nearest neighbors is found to carry identical topological charges. The resulting short- to mid-range 
order among the topological charges is reminiscent of Debye screening clouds in a plasma. These ordering patterns 
can be traced to the dependence of the underlying interactions on the relative color orientations between neighboring 
dimerons. Ensemble and configuration averages of these color orientations indicate, on the other hand, that in 
the mean the intermediate-range attraction between neighboring dimerons and anti-dimerons is almost completely 
compensated by their shorter- and longer-range repulsion and by the ensemble entropy. In any case, with increasing 
coupling and hence entropy the observed topological order weakens rather rapidly, together with the underlying color 
correlations. 

The most far-reaching impact of the increasing bare coupling, however, manifests itself in the structural changes 
which it induces in the individual dimerons. In particular, we found the first robust evidence, established in fully 
interacting ensembles, for the dimerons to gradually dissociate into their meron constituents when the coupling grows. 
This process is driven by a competition between the coupling-dependent interactions and the entropy. Several com- 
plementary pieces of evidence further indicate that the meron constituents continuously replace the dimerons as the 
dynamically relevant degrees of freedom. With growing coupling the interactions between merons belonging to dif- 
ferent dimerons therefore increasingly determine the ensemble properties. This renders the dimerous' meron partners 
more and more independent (except for their frozen relative color orientation) and allows our field configurations to 
approximate more closely those of meron ensembles. In particular, however, the changing nature of the dynamically 
active topological-charge carriers and of their coupling-dependent correlations leaves various revealing imprints on 
bulk ensemble properties. 

Among the latter, the topological susceptibility is particularly relevant in our context. When calculated in a 
sequence of growing evaluation volumes up to the largest accessible ones, it is found to level off towards the end of the 
sequence. This indicates or at least foreshadows the expected saturation of the topological-charge fluctuations in the 
thermodynamic limit. Moreover, the resulting values of the susceptibility in the intermediate, physical coupling range 
remain within errors practically constant. This behavior may be related to the similarly weak temperature dependence 
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of the topological susceptibility in Yang-Mills theory below the deconfinement temperature. Our intermediate-eoupling 
results reproduce those of comparable single-meron ensembles rather closely, furthermore, which seems to indicate 
that (on average relatively loosely bound) meron centers of dimerous and single merons produce topological charge 
fluctuations of roughly equal strength. On the other hand, our results for the topological susceptibility underestimate 
lattice and instanton ensemble values by about 50%, perhaps due to insufficient saturation in our evaluation volumes 
and the underlying boundary effects. In any case, the dimcron ensemble predictions for the susceptibility could be 
straightforwardly adjusted by increasing the dimeron density. 

Since the long-range color tails of less bound meron centers disorder the vacuum more strongly, we have payed 
particular attention to the confinement properties of the dimeron ensembles. The latter were monitored by evaluating 
ensemble averages of rectangular Wilson loops and by extracting the associated static quark-antiquark potentials for 
a wide range of coupling values. Initially, the resulting potentials are found to rise quadratically with the interquark 
separation. This behavior is typical for (sufficiently dilute) instanton ensembles and indicates that our short-distance 
potentials are dominated by contributions from moderately dissociated dimerous which act essentially as singular- 
gauge instantons. For quark-antiquark separations beyond about half a fermi, on the other hand, the Wilson loops 
indeed develop an area-law behavior and the potentials consequently show the linear distance dependence characteristic 
of quark confinement. The associated string tension is found to grow with the coupling and to reach values of the 
experimentally expected order of the magnitude in the physical coupling range. 

Hence we have shown - without reliance on weak-coupling or low-density approximations - that sufficiently strongly 
coupled and therefore on average far enough dissociated dimeron ensembles indeed produce confining long-range 
correlations of about the physically required strength. Since the latter originate from the meron tails of the more 
strongly dissociated dimerous, these results also provide complementary evidence for single-meron ensembles to confine 
(even though our finite simulation volume and our limited field configuration space prevent the dimerons from breaking 
up completely). Moreover, both dimeron and single-meron ensembles are found to generate confining static-quark 
potentials of comparable strengths at our largest accessible interquark distances. Since the linear rise of the meron- 
induced potential turns out to continue beyond these distances, furthermore, our above reasoning indicates that the 
same should happen in dimeron ensembles. Hence we expect an ongoing linear rise of the dimeron-induced potentials 
into the distance region where the underlying color flux tubes are fully developed. 

While dimeron ensembles thus share with single-meron ensembles the main features of the confinement mechanism, 
they additionally provide a more complete description of the Yang-Mills vacuum. The first evidence is that our 
results establish how a confining meron component in the vacuum can robustly emerge from instanton and dimeron 
dissociation. This provides, in particular, the first quantitative confirmation for a key ingredient of the meron-induced 
confinement scenario envisioned more than three decades ago on the basis of qualitative arguments. In addition, 
our results on the coupling dependence of dimeron ensemble properties should contain useful information on their 
temperature dependence and, in particular, on the impact of temperature-induced variations of the gauge coupling 
(although we did not pursue this issue in the present paper) . 

Furthermore, we found evidence for the physical dimeron density to be about twice as large as typical instanton 
densities. This allows confining dimeron ensembles to retain a phenomenologically relevant fraction of contracted and 
thus more weakly interacting, "instanton-like" dimerons. An independent indication for their considerable density is 
the strong quadratic rise of our heavy-quark potentials at short interquark distances. The instanton-like component 
should adapt to the presence of light quarks in essentially the same fashion as in phenomenologically successful 
instanton vacuum models, furthermore, since the underlying quark zero-modes in the instanton background are of 
topological origin and therefore little affected by the moderate deformations of instantons into contracted dimerons. 
Hence, in addition to confining color, dimeron ensembles should be able to reproduce much of the successful vacuum 
and hadron phenomenology predicted by (non-confining) singular-gauge instanton ensembles. 
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[68] As in several lower-dimensional systems, in the thermodynamic limit the dimeron dissociation process may therefore lead 

to a sharp transition of Kosterlitz-Thouless type [s^, [s^l from a deconfining to a confining phase. 
[69] i.e. the 2nd Chern class of the associated principal fiber bundle 

[70] More generally, these are the N'^ — 1 "Euler angles" which parametrize (constant) SU(A'^c) group elements. 
[71] It may be useful to recall that the superposition of singular-gauge dimerons is not gauge-equivalent to a superposition of 
regular-gauge dimerons. 

[72] From the practical point of view this is fortunate since no multi-meron solutions with |Q| > 1 seem to be known analytically. 
[73] More precisely, all four coefficients of the quaternion representation are varied independently. The constraint among 

them is reinstated only afterwards by rescaling all coefficients with a common factor. (Ifence (ci) represents the Haar 

measure on SU(2) ~ 5"^.) 

[74] The singular-gauge induced singularities are so strongly localized (even after regularization, cf. Sec. Ill Cp . on the other 
hand, that they have practically no impact on the simulations. 

[75] In practice, these densities were calculated by dividing the box volume into 60 concentric, spherical shells around the 
reference pseudoparticle, counting the number of either equally or oppositely charged pseudoparticle centers contained in 
it, and dividing by the shell volume Vshcii. The latter is calculated analytically as long as the shells lie fully inside the 
ensemble volume (cf. Sec. lIII A| . and otherwise as Vshcii = VNahsu/Nv where Nv points are randomly distributed over the 
ensemble box and A'shoii is the number of them which fall into the shell. 

[76] This screening behavior should not be confused with the light-quark and anomaly-induced topological charge screening in 
the QCD vacuum [43], with its strong impact on rj' meson 48] and pseudoscalar glueball properties 49, 5Cf. 

[77] In order to reduce the impact of boundary effects we include only those (anti-)dimerons whose centers lie inside a (hy- 
per)sphere of radius one around the mid point of the box. The inter-meron separation 2 (|a|) of the dimerons then restricts 
the maximal distance of the meron center from the box center to 1 -f (|a|), i.e. on average meron centers lie maximally a 
distance (|a|) outside the ball. For = 1, where the boundary efi'ects are strongest, the smallest (|a|) ~ 0.27 then reduces 
their impact by ensuring that the meron centers remain closest to the sphere. For the maximal couplings, on the other 
hand, one finds (|a|) ~ 3/2. 

[78] Even for Nd = Njy strict CP invariance would hold only in the infinite- volume limit and/or with infinite ensemble statistics, 
incidentally, since the topological charge density of the dimerons is not fully contained in our evaluation volumes. However, 
this boundary artefact is mitigated already in relatively small volumes by the short-range anticorrelations among the 
dimerons' topological charges (cf. Sec. IIVB|) . 

[79] Our V resolution is obtained by using volumes with extent 1.2 and step size 0.1 in the three spacial directions, and a step 
width of 0.03 in the temporal direction which maximally extends to 4.0. 

[80] The values of very small Wilson loops could alternatively be obtained analytically since our gauge field configurations (|13|) 
become constant at length scales <g p (cf. Ref. [3)- 

[81] The statistical error of the data points from the subvolumina is somewhat smaller than that from their central hyperplanes. 
When plotting {W [C {R, 2R)]) as a function of R, incidentally, some seemingly correlated discrepancies between these data 
sets become nevertheless visible for = 1 but not for larger couplings. Hence these distortions may indeed originate from 
boundary effects (cf. Sec. IIIIDI) and appear to be sufficiently suppressed in the chosen subvolumina. 

[82] Linear short-distance contributions to the static quark potential (originating e.g. from a nonlocal two-dimensional gluon 
condensate and modeled by a tachyonic UV gluon mass) were argued for as well (It is conceivable that the latter 
may also contribute in dimeron ensembles where the linear behavior of V (R) sets in at rather small distances.) 

[83] In Refs. [l3, [3 argued that the meron peak contributions, which diverge logarithmically for p — )■ 0, should be 

subtracted from (s) in meron ensembles. The remaining, about 25-35% reduced "background" action density is then 
considered as the physical result, renormalized at . For our more localized singular-gauge dimerons this procedure will 
be less effective, and we shall not attempt to implement it here. 

[84] The relation given in Ref. contains an error. We thank Thomas Schafer for clarifying correspondence on this issue. 

[85] Even at finite coupling, correlations between instantons are known to be fairly weak in pure Yang-Mills theory and in 
singular-gauge instanton ensembles (without light quarks). Hence the topological susceptibility can be used to estimate 
the instanton density in the Yang-Mills vacuum [2^ . 

[86] Indications for the coexistence of instantons and merons at comparable densities in the Yang-Mills vacuum were also found 
in a variational approach i62|] . 
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FIG. 1: The multi-layered multigrid designed to control boundary effects when numerically evaluating the action density. 
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FIG. 2: The thermalization histories of simulations with — 1 (left panels) and — 25 (right panels) as a function of the 
number nswp of sweeps. The uppermost row shows the evolution of the maximal step widths in position (lower, black curves) 
and color (uper, red curves) space. In the second row we plot the average distances d between nearest-neighbor dimerous of 
equal (black curves, uppermost for = 1), opposite (red curves, intermediate for g^ = 1) and arbitrary (lowest-lying, blue 
curves) topological charge. The last row shows the average probability / for the nearest neighbor of a given meron center to 
have opposite topological charge. (The dotted vertical lines indicate the sweep numbers at which the resolution of the action 
sampling is increased. For a stringent test of the resolution quality we enhance the corresponding action changes by smearing 
the meron-center singularities less broadly, corresponding to — 0.2, in the = 1 simulation (left panels). For g^ = 25 (right 
panels) we use our standard value = 0.3.) 
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FIG. 3: The averaged action density (s) (black squares, full line) and quadratic "plaquette" Wilson loop expectation values 
(W) (red circles, dotted line, multiplied by 450) along the time direction (left panel) and along the spacial diagonals (right 
panel) from the center of the simulation box for = 1. (The dashed horizontal lines are fits to constant plateaus for t,r < 1. 
The two vertical lines indicate the boundaries of the core and sampling volumes.) 
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FIG. 5: The average separation 2 {|a|) between the (anti-) meron centers of the (anti-) dimerons as a function of the squared 
gauge coupling. (The error bars are smaller than the plot symbols and the dotted vertical line indicates the scale change at 
= 1000.) 
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FIG. 6: The action (left panels) and selfduality (right panels) densities of a typical = 1 ensemble configuration in two cross 
sections parallel to the x\ — X4 plane. Lighter shades of gray indicate larger values. The ordinates denote the xi and the abcissas 
the Xi direction. The upper row shows the densities in the hyperplane at {x2,X3) = (0.4,0.4) and the lower row those in the 
plane at {x2,X3) = (0.4, 1.2). 
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FIG. 7: The average radial density of topologically equally (black squares, full line) and oppositely (red bullets, dotted line) 
charged dimerous pj-, (left panel) and merons pj^j (right panel) as a function of the distance from a fixed (anti-) (di)meron at 
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FIG. 8: The same as in Fig. [71 but for = 25. 
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FIG. 9: The same as in Fig. but for = 100. 
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FIG. 10: The average probability (/□□) for the nearest neighbor of a dimeron to have opposite topological charge. (The error 
bars are smaller than the plot symbols and the dotted vertical line indicates the scale change at — 1000.) 




FIG. 11: The left panel shows the average distance (d) of a reference dimeron center from its nearest neighbor with equal (black 
squares, full line) and opposite (red bullets, dotted line) topological charge Q. The nearest-neighbor distance (d) independent 
of Q is also included (blue crosses, dashed line). The right panel shows the same curves for the average distances between 
meron instead of dimeron centers. (The error bars are smaller than the plot symbols and the dotted vertical line indicates the 
scale change at = IQQQ.) 
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FIG. 12: The distribution of the average color angle (q) between nearest neighbor pseudoparticles of equal (black squares, 
full hne) and opposite (red bullets, dotted line) topological charge at — 1. For comparison, we also show (a) for a random 
pseudoparticle distribution (i.e. at = oo; blue crosses, dashed line) which is independent of the relative topological charge 
of the neighbors. (The error bars indicate the statistical standard deviation of the averages per bin.) 
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FIG. 13: The standard deviation of the (q) distribution between nearest neighbors of equal (black squares, full line) and 
opposite (red bullets, dotted line) topological charge Q as a function of g^. The standard deviation of the color orientation for 
nearest neighbors independent of Q is also included (blue crosses, dashed line). 
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FIG. 15: The logarithm of the expectation value of rectangular Wilson loops with area A (and side length 7? = T/2 = 
for the square coupling values — 1, 25, 100, 1000 and oo. (For increasing coupling the data lie below each other.) The black 
squares (boxes) are obtained from the loops in the full evaluation volume, the red dots (asterisks) from the loops in the central 
hyperplanes. The fit curves to / (A), as described in the text, are drawn as full lines in their fit intervals. 
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FIG. 16: The heavy-quark potentials V{R) for the square coupling values — 1, 25 and 100. (The stronger potentials 
correspond to larger couphngs.) 
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FIG. 17: The dimensionless ratio cr^^'^/x^^'^ ^ function of the square coupling g^. 



